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The  Co-60  gamma  radiolysis  of  gaseous  propane  was 
studied  at  100  torr  pressure  and  25  °C,  both  pure  and  with 
10%  added  oxygen.  In  the  unscavenged  system  the  major 
products  and  their  G-values  were  hydrogen,  4.99;  methane, 
1.30;  ethane,  1.95;  iso-butane,  0.61;  n-butane,  0.25;  i- 
pentane,  0.42;  n-pentane,  0.14;  and  hexanes,  0,89.  Minor 
products  were  heptanes,  0,082;  octanes,  0,067;  nonanes, 
0.088,  and  decanes,  0.033.  Small  yields  of  ethylene  and 
propylene  were  also  observed.  Yields  in  the  presence  of 
10%  added  oxygen  were  hydrogen,  1.87;  methane,  0.83;  and 
ethane,  1.22.  Higher  saturated  hydrocarbons  were 
eliminated. 

In  the  presence  of  10%  oxygen,  the  major  oxygen  con- 
taining products  and  their  G-values  were  acetone,  0.98; 
iso-propyl   alcohol,    0.86;      propionaldehyde,    0.43;  n-propyl 
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alcohol,  0.11;  acraldehyde,  0.14;  and  allyl  alcohol,  0.038. 
Minor  products  were  iso-butyl  alcohol,  0.0085;  n-butyl 
alcohol,  0.0029;  tert-amyl  alcohol,  0.0074;  iso-amyl  alco- 
hol, 0.0013;  and  n-amyl  alcohol,  0.0016.  When  the  oxygen 
concentration  was  decreased  below  5%,  it  was  found  that  the 
yields  of  acetone,  iso-propyl  alcohol  and  n-propyl  alcohol 
increased,  propionaldehyde  yield  decreased  and  yields  of 
other  products  remained  constant. 

The  reaction  scheme  for  formation  of  major  products 
was  examined  using  computer  modeling  based  on  24  reactions 
in  the  unscavenged  system  and  28  reactions  in  the  propane- 
oxygen  system.  Yields  could  be  brought  into  agreement  with 
the  data  within  experimental  error  in  nearly  all  cases,  but 
in  the  pure  propane  system  it  was  necessary  to  assume  that 
the  molecular  hydrogen  yield  was  accompanied  by  the  deposi- 
tion of  polymer  on  the  vessel  wall. 
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I.  INTRODUCTION 


Production  of  useful  compounds  by  the  chemical  reac- 
tions taking  place  under  the  action  of  radiation  is  one  of 
the  important  branches  of  atomic  energy  application.  The 
radiation  chemistry  of  propane  is  of  special  importance 
because  it  is  the  lowest  normal  alkane  which  provides  both 
primary  and  secondary  carbon  atom  sites.  There  have  been  a 
number  of  research  papers  published  concerning  the 
radiolysis  of  gaseous  propane  both  with  and  without 
scavengers  (1-23). 

In  1926,  Lind  and  Bardwell  (1)  reported  the  formation 
of  product  as  heavy  as  pentane  during  the  bombardment  of 
propane  gas  with  6-Mev  alpha  particles  from  radon  gas. 

In  1952,  Gevantman  and  Williams  (2)  made  a  study  of 
the  free  radicals  and  their  relative  production  in  the 
radiolysis  of  propane  with  x-rays  and  electrons,  but  they 
did  not  report  a  complete  product  distribution.  According 
to  them  in  the  radiolysis  of  propane,  there  was  a  distinct 
effect  of  the  propane  pressure.  With  increasing  pressure 
the  relative  yield  of  the  methyl  decreased  and  the  yield  of 
methylene  increased  in  such  a  fashion  that  the  sum  of  the 
yields  of  methyl  and  methylene  was  approximately  constant, 
and  the  yields  of  other  radicals  remained  approximately 
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constant.  To  explain  these  two  proposals  were  put 
forward: 

1)  Reaction  of  thermal  methyl  radicals  with  propane 
might  produce  methylene  radicals. 

CH3*  +  C3H8  — >  CH2  +  H2  +  C3H7* 

Since  such  reactions  must  be  highly  endothermic,  this 
was  not  a  satisfactory  proposal. 

2)  Methyl  and  methylene  radicals  are  both  primary 
products  of  excited  propane 

C3H8*  — >  CH3*  +  C2H4  +  H* 
— >  CH2    +  C2H6 

A  pressure  effect  on  the  relative  efficiencies  of 
these  two  processes  was  assumed. 

In  1959,  Yang  and  Manno  (4)  carried  out  the  radiolysis 
of  propane  using  gamma  radiation  and  reported  the  identi- 
fication of  ten  different  compounds  ranging  from  hydrogen 
to  pentanes.  After  irradiation,  the  products  were  analyzed 
by  gas  chromatography.  They  also  studied  the  effect  of  NO 
as  a  scavenger  on  the  yields  of  various  products.  Assuming 
the  G(H2)  in  propane  to  be  5.9,  the  product  distribution 
relative  to  hydrogen  yield  was  reported. 

In  1960,  Birdwell  and  Crawford  (5)  determined  the 
primary  products  and  their  relative  distribution  in  the 
gamma   radiolysis  of  propane  gas  at  low  irradiation  doses. 
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They  identified  compounds  including  hydrogen,  hexanes, 
ethane,  methane,  iso-butane,  isopentane,  n-butane,  and  n- 
pentane;  alkenes,  such  as  propylene,  ethylene  etc.  were 
indicated  in  the  mass  spectra  of  the  irradiated  propane 
samples,  but  their  yields  were  too  small  for  identification 
by  gas  chromatograph.  These  experiments  were  carried  out 
at  a  pressure  of  75  psi  (388  cm)  at  a  temperature  of 
approximately  25  °C  using  a  Co-60  gamma  source. 

Birdwell  and  Crawford  (5)  were  the  first  to  identify 
hexanes  in  the  gamma  radiolysis  of  propane.  They  reported 
the  product  distribution  relative  to  the  hydrogen  yield. 
No  experiments  were  done  by  them  using  a  scavenger.  No 
mechanism  was  given  by  Birdwell  and  Crawford  or  by  any 
other  previously  mentioned  authors  for  the  formation  of  the 
various  products  in  the  radiolysis  of  gaseous  propane. 

In  1960,  Back  (6)  carried  out  the  gamma  radiolysis  of 
hydrocarbon  gases  at  very  low  conversions  and  used 
ethylene/propylene  as  scavengers  for  hydrogen  atoms  and 
concluded  that  two-thirds  of  the  hydrogen  produced  in  the 
gamma  radiolysis  of  ethane,  propane,  n-butane  and  n-pentane 
in  the  gas  phase  came  from  hydrogen  atoms. 

Among  the  most  recent  publications,  the  papers  of 
Fujisaki  et  al.  (22)  and  Gawlowski  et  al.  (23)  deal  with 
special  topics  of  primary  processes  of  the  hydrogen  yield 
and  impurity  effects,  respectively. 
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Many  other  previous  workers  on  the  propane  system  had 
provided  information  on  specific  aspects  of  the  chemistry 
involved,  including  ionic  fragmentation  of  propane  (24-28); 
ion-molecule  reaction  chemistry  (29-33);  and  decomposition 
of  excited  propane,  studied  by  VUV  photochemistry  (14,34). 
Numerous  ionic  processes  occur  along  with  the 
radical/excited  molecular  reactions  which  contribute 
substantially  to  certain  product  yields.  A  sequence  of 
elementary  reactions  was  observed  which  lead  to  C3Hg"*'  and 
C3H^'''  as  stable  ions  (33).  In  the  mass  spectrum  of  propane 
three  ions  constitute  over  60%  of  the  total  ionization  (18) 
C2H5+  (30.5%),  C2H4+  (13.1%)  and  C2H3'*"  (12.3%).  The 
proposed  route  is 

^3^8^  ~>  C2H5"''  +  CH3' 
~>  C2H4+  +  CH4 
~>  C2H3'*'  +  H2  +  CHj* 

The  work  of  Bone,  Sieck,  and  Futrell  (20)  gives  a 
broad  overview  of  propane  radiolysis  with  emphasis  on  ionic 
processes  and  noble  gas  sensitization. 

In  the  analysis  of  oxygenated  products  it  is  necessary 
to  find  out  what  products  are  formed  in  the  pure  propane 
system  and  how  the  yields  of  these  products  are  affected 
when  a  scavenger  such  as  oxygen  is  added.  It  would  be 
useful  to  report  our  data  and  compare  it  to  previous  work 
and  carry  out  a  computer  analysis  of  the  proposed  reaction 
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scheme,  which  may  be  of  interest  in  relation  both  to  the 
earlier  work  and  to  our  subsequent  investigation  of  the 
oxygen  scavenged  system  (35). 

Oxygen  has  frequently  been  used  as  a  free  radical 
scavenger  in  the  gamma  radiolysis  of  both  liquid  and  gas 
phase  hydrocarbon  systems  (36).  Most  workers  have  recorded 
only  the  effect  of  oxygen  in  decreasing  the  yields  of 
various  hydrocarbon  products  formed  in  the  unscavenged 
system.  It  is  generally  assumed  that  the  portion  of  each 
product  yield  eliminated  in  this  way  is  due  to  free  radical 
reactions.  In  most  instances,  no  attempt  was  made  to 
characterize  and  measure  the  oxygen-containing  products 
which  are  formed  in  the  presence  of  oxygen.  Molecular 
oxygen  is  an  efficient  acceptor  for  most  organic  radicals. 
A  study  of  radiation  chemical  processes  in  organic  systems 
in  the  presence  of  molecular  oxygen  allows  one  to  ascertain 
the  role  of  primary  radicals  in  the  basic  mechanism  of 
these  processes,  and  at  the  same  time  opens  new  ways  of 
producing  oxidation  products  of  practical  interest  (37). 
Some  work  has  been  published  on  the  radiolytic  oxidation  of 
liquid  hydrocarbons  (38-51)  as  well  as  on  gases  (51-60), 
especially  methane. 

In  1926,  Lind  and  Bardwell  (1)  irradiated  mixtures  of 
methane  and  oxygen  with  radon  alpha  rays.  The  end  products 
of   the   reaction  were   reported   to  be  carbon  dioxide  and 
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water.  No  attempt  was  made  to  identify  any  intermediate 
products,  nor  to  establish  a  reaction  mechanism. 

In  1958,  Mikhailov,  Kyimora  and  Bogdonov  (38)  studied 
the  irradiation  of  methane-oxygen  mixtures  with  a  beam  of 
electrons  produced  at  112  kV,  The  oxygenated  products 
identified  were  carbon  monoxide,  carbon  dioxide,  and  water. 
The  presence  of  peroxide,  aldehydes,  acids,  and  alcohols 
was  shown  by  group  reaction  but  these  products  were  not 
characterized. 

In  1959  Liebenthal,  Albright,  and  Sesonske  (43)  irra- 
diated mixtures  of  propane  and  oxygen  in  the  temperature 
range  of  350  °C  to  395  °C  at  atmospheric  pressure.  The 
propane-to-oxygen  mole  ratios  were  3:1,  5:1,  and  10:1.  The 
yields  of  formaldehyde,  acetaldehyde  and  methanol  increased 
significantly  under  gamma  irradiation,  sometimes  as  much  as 
100%   -  200%. 

In  1959  Bakh  and  T'ung  T'ien-chen  (52)  investigated 
the  formation  of  liquid  products  in  the  oxidation  of  n- 
heptane  and  n-nonane  using  x-rays  over  a  temperature  range 
of  -80°C  to  +130°C.  The  acids  formed  were  estimated  by 
titration  with  sodium  methylates  in  a  nonaqueous  medium. 
The  scheme  for  the  production  of  acid  was  as  follows: 


H 


RCHO2  > 


RCOOH 


R 
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In  1961,  Johnson  and  Salmon  (39)  reported  the  forma- 
tion of  methyl  hydroperoxide  in  the  radiolysis  of  methane 
in  the  presence  of  oxygen. 

In  1962,  Slavinskaya,  Zhitneva,  Kamens tskaya ,  and 
Pshezhetskii  (46)  described  the  kinetics  of  the  chain  oxi- 
dation of  butane  with  oxygen  under  the  action  of  fast 
electrons,  and  the  characteristic  features  of  the  reaction 
mechanism,  particularly  the  effect  of  the  radiation  on  the 
branching  of  the  reaction  chain.  The  peroxides  and  alde- 
hydes were  determined  polarographically.  The  scheme  is 
based  on  processes  which  are  most  probable  in  thermal 
oxidation. 

RH  — >  R*  +  H* 

R*  +        — >  ROj* 

R02*   +  RH  — >  ROOH  +  R* 

RO2*   — >  R^HO  +  RjO* 

In  1965,  Knox  (47)  proposed  a  mechanism  for  the  low 
temperature  oxidation  of  alkanes  at  300  to  400  °C  where 
HO2.  is  converted  to  OH.  radicals.  This  conversion  involves 
a  co-oxidation  of  the  olefin  AB,  which  is  the  primary 
oxidation  product  of  most  alkanes.  The  key  to  the  conver- 
sion mechanism  is  the  sequence  shown: 

HO2                 O2  HO2 
AB   >  ABOOH   >  OOABOOH   >  HOOABOOH   >   20H  +  AO  +  BO 
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In  1966,  Hummel  and  Hearne  (40)  studied  the  radiolytic 
oxidation  of  methane.  They  reported  the  initial  G-values  of 
thirteen  oxygenated  and  four  hydrocarbon  products  for  a 
mixture  of  10%  O2  in  methane  over  a  range  of  pressures  (1 
to  76  cm  Hg)  and  at  two  dose  rates. 

In  1968,  Sochet,  Sawerysyn,  and  Lucquin  (48)  studied 
the  reaction  of  propane  and  oxygen  at  various  pressures, 
temperatures,  and  concentrations.  The  analysis  of  some 
reaction  products  indicates  an  increase  in  the  amounts  of 
methane,  ethane,  acetaldehyde,  ethyl  alcohol,  propyl  alco- 
hol, and  especially  iso-propyl  alcohol  and  a  decrease  in 
the  formation  of  hydrogen  peroxide  and  olefin.  All  these 
results  can  be  explained  by  radical  reactions  such  as 

R*    +   R02*    (or   R02*)  >    ROOR  >  2R0* 

 >  oxygenated  products  and  R*  +  R*  >  rr 

In  1978,  Kunjappu  and  Rao  (60)  studied  the  radiation 
induced  oxidation  of  n-heptane  at  30  °C  using  a  Co-60  gamma 
source.  They  measured  the  yields  of  heptanols  and  hepta- 
nones;  no  acids  were  observed.  Since  4-,  3-,  and  2-hepta- 
nones  are  the  only  carbonyl  compounds  produced  during  the 
oxidation  of  n-heptane,  this  result  suggests  that  added 
oxygen  interferes  with  reaction  pathways  which  lead  to  C-C 
bond  rupture  in  an  unscavenged  system.  Carbonyl  and  alco- 
holic products  are  formed  by  the  decomposition  of  peroxy 
radicals     through  a  disproportionation  reaction  between 
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them.  Iodine  and  triphenyl  methane  (TRITAN)  were  used  as 
scavengers.  The  scavenging  reactions  were  represented  as 
follows : 

^7^16        ^  ^7^15*  ^' 

^7^15*  +  02   >  C7H1502- 

07^^1502*  +  ^1^15^2'  >  CyHj^jOH  +  C7H2^40  +  O2 

07^15*       Scavenger  >  products 

According  to  Kunjappu  and  Rao,  a  fraction  of  the 
heptanones  is  formed  by  the  direct  interaction  of  the 
excited  state  with  molecular  oxygen. 

.         *  °2 

 >  C7H24O  +  H2O 

Below  is  a  brief  outline  of  the  modern  concept  of  the 
mechanism  of  hydrocarbon  oxidation  in  the  gas  phase 
(36,61);  the  reaction  pre-supposes  the  presence  of  alkyl 
radicals  which  are  usually  formed  via  thermal  decomposition 
of  an  initiator.  Oxygen  readily  adds  to  free  radicals  and 
if  present  will  almost  invariably  play  a  part  in  radiation 
induced  reactions.  The  combination  of  oxygen  with  radicals 
is  easily  understood  if  molecular  oxygen  is  written  as  a 
diradical  'O-O*  in  recognition  of  its  triplet-ground  state 
(36).     The  product  is  a  relatively  stable  peroxy  radical. 

R*  +  O2  — >  R-O-0*  (a) 
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A  peroxy  radical  may  abstract  hydrogen  from  the  substrate 

R02*  +  RH  — >  RO2H  +  R*  (b) 

Reactions  a  and  b  together  constitute  a  chain  reaction 
of  the  type  responsible  for  the  autoxidation  of  a  variety 
of  organic  compounds,  Autoxidation  is  the  slow  oxidation 
of  an  organic  compound  by  oxygen. 

Combination  of  two  peroxy  radicals  may  give  an  or- 
ganic peroxide 

2R02*  — >  R-O-O-R  +  O2 
but  the  products  more  often  are  two  alkoxy  radicals. 

2R02*  — >  2R0*  +  O2 

At  higher  temperatures 

R02*   +  RH   >  ROOH  — >  RO*   +  OH* 

Alkoxy  radicals  are  more  reactive  than  peroxy  radicals 
and  can  abstract  hydrogen  from  most  organic  compounds 

RO  +  RH  — >  ROH  +  R* 

The  combination  of  two  alkoxy  radicals  to  give  a 
peroxide  appears  improbable  relative  to  disproportionation 
to  give  an  alcohol  and  carbonyl  compound 

2RCH20*   — >  RCH2OH  +  RCHO 
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similarly,  reactions  between  alkoxy  radicals  and  an 
alkyl  radical  are  less  likely  to  give  an  ether  (RCHjOR) 
than  a  carbonyl  compound  and  a  hydrocarbon. 

RCH20*   +  R*   >  RCHO  +  RH 

According  to  Benson  (62),  olefins  turn  out  to  be  more 
easily  oxidized  than  alkanes.  Allylic  H  atom  abstraction 
leads  to 

CH3CH=CH2  +  R0«   >  ROH  +  CH2=CHCH2 

CH2=CHCH2  +  O2   >  CH2=CHCH200* 

We  have  selected  propane  as  the  most  appropriate  sub- 
strate for  an  investigation  of  radiolytic  oxidation  because 
it  is  the  simplest  species  providing  both  primary  and 
secondary  carbon  atom  sites.  Since  no  work  was  done  pre- 
viously to  identify  the  products  in  propane-oxygen  system, 
we  thought  it  would  be  useful  to  identify  the  gamma  ray 
oxidation  products  of  gaseous  propane  and  carry  out  a 
computer  analysis  of  the  proposed  reaction  scheme. 
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II.     EXPERIMENTAL  METHODS  AND  EQUIPMENT 

A.  Reagents  and  Their  Purifications 

Propane;  Research  grade  (99.99%)  propane  obtained  from 
Matheson  Company  was  used.  It  was  degassed  by  several 
f reeze-pump-thaw  cycles  and  stored  in  a  vessel  attached  to 
the  vacuum  line.  This  material  showed  no  observable 
impurities  by  flame  ionization  gas-chromatography. 

Ethylene;  Matheson  Company  CP.  grade  (99.5%)  ethylene 
was  passed  through  a  BaO  drying  tube  into  a  storage  vessel 
on  the  vacuum  line.  Possible  air  contamination  was  then 
removed  by  several  f reeze-pump- thaw  cycles.  Ethylene 
showed  no  observable  impurities  by  flame  ionization  gas- 
chromatography. 

Oxygen;  Matheson  Company  Extra  Dry  (99.6%)  oxygen  was 
used  as  received  from  the  manufacturer. 

Other  Reagents;  Chemicals  used  for  chromatographic 
calibration  standards  were  used  as  received. 
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B.     Sample  Preparation 

Vacuum  System;     All  samples  for  radiolysis  were  pre- 
pared on    a  vacuum  line.    The  vacuum  system  used  is  shown 
in  Figures  1  and  2.  The  pumping  system  was  a  Welch  Duo-Seal 
pump  connected  through  a  liquid  nitrogen  trap  to  a  two 
stage  mercury  diffusion  pump.     These  pumps  were  connected 
to  the  main  manifold  through  a  second  liquid  nitrogen  trap 
(X)  and  a  stopcock   (V^^).     Attached  to  the     main  manifold 
were  the  following:     a  mercury  manometer   (M)  ,   a  vacuum 
thermocouple  gauge   (G-^) ,    inlets  for  propane  and  oxygen  (I^ 
and        ,  two  storage  bulbs  (B^  and  Bj)  and  a  Toepler  pump- 
McLeod  gauge  apparatus  Fig.  2.    Teflon  stopcocks  (T)  were 
used    in    the    submanifold.      Two    stopcocks    isolated  the 
Toepler  pump  from  the  main  manifold.     A  submanifold  was 
provided  to  attach  the  radiolysis  vessel   (R)   and  sample 
loop  (Sj)  to  the  input  of  the  combined    Toepler  pump-McLeod 
gauge  system,    via  a  liquid  nitrogen  U-trap  packed  with 
glass  helices;   this  enabled  the  transfer  of  the  radiolysis 
products  to  the  sample  loop  (Sj^)   for  analysis  following 
irradiation.     A  special   three-way  stopcock  allowed  the 
Toepler   pump   to    transfer   non-condens i bl e  radiolysis 
products  either  into  a  series  of  calibrated  volume  (thus 
acting  as  a  McLeod  gauge)   or  into  a  sample  loop  (S2)  for 
analysis   by   gas-chromatography;    an   additional,  custom 
tabulation     on  the  base  of  the  vacuum  stopcock  allowed 
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evacuation  of  the  sample  loop  (S2).  This  system  was  also 
equipped  exclusively  with  Teflon  stopcocks  (T)  and  had  a 
separate  thermocouple  vacuum  gauge  (G2)  for  monitoring 
sample  transfer  operations. 

Propane;  Before  each  sample  preparation,  propane  was 
transferred  from  cylinder  (C-j^)  to  the  storage  reservoir 
(B]^)  and  was  de-aerated  by  repeated  f reeze-pump-thaw  cycles 
until  the  thermocouple  gauge  indicated  the  absence  of  air. 
Valve  T3  was  then  closed  and  the  condensed  gas  was  brought 
to  room  temperature,  Radiolysis  vessels,  having  known 
volumes,  were  attached  to  0^  and  Valve  T^  was  closed 

and  valve  T3  was  opened.  The  sample  was  allowed  to  expand 
into  the  radiolysis  vessels  (R).  During  the  course  of  this 
operation,  the  pressure  was  monitored  with  a  mercury 
manometer  (M),  At  the  desired  pressure,  valve  (T3)  and  the 
valves  of  the  radiolysis  vessels  were  closed  and  the  excess 
material  was  condensed  back  into  the  storage  vessel  with 
liquid  nitrogen  by  opening  valve  T3,  Knowing  the  tempera- 
ture, pressure,  and  volume,  the  ideal  gas  law  was  used  to 
calculate  the  amount  of  sample  in  the  radiolysis  vessels 
(R)  . 

Propane  with  added  Oxygen;  Propane  gas  stored  in  the 
storage  vessel  (B-j^)  was  de-aerated  by  repeated  freeze-pump- 
thaw  cycles  and  was  brought  to  room  temperature,  A  known 
amount  of  oxygen  was  taken  in       and  then  propane  was  added 
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to  ^^'^  sample  was  allowed   to   mix.    The  desired 

amount  of  mixture  was  allowed  to  expand  into  the  previously 
evacuated  radiolysis  vessels  which  were  attached  to  0-^  and 
To  ensure  that  the  right  amount  of  oxygen  was  taken, 
the  reaction  vessel  was  attached  to  the  sample  loop  which 
was  attached  to  the  Toepler  pump-McLeod  gauge  system  as 
shown  in  Fig.  2.  The  cold  trap  and  the  sample  loop  were 
immersed  in  liquid  nitrogen  .  Valve  7^,2  was  closed  and  the 
valve  of  the  reaction  vessel  was  opened.  Propane  was 
allowed  to  condense  in  the  sample  loop  and  cold  trap.  The 
amount  of  oxygen  was  determined  using  the  Toepler-McLeod 
apparatus.  The  amount  of  oxygen  was  also  determined  by  gas 
chromatography  after  transferring  the  oxygen  to  the  sample 
loop  Sj. 

C.     Sample  Irradiation 

Radiation  Source  and  Vessel;  Irradiations  were  car- 
ried out  using  a  300  curie  Co-60  gamma  ray  source,  which 
has  been  described  previously  (63).  Figure  3  is  a  cross- 
sectional  view  of  the  irradiator.  The  decay  scheme  of  Co-60 
is  shown  in  Fig.  4. 

The  helium  arc  welded  nickel  radiolysis  vessel  used 
is  shown  in  Fig.  5.  Two  vessels  were  used,  having  volumes 
108.5  cc.  and  112.0  cc,  respectively.  The  vessels  were 
made  of  pure  nickel  having  Hoke  monel  diaphragm  valves 


Legend:  (A)  counterweight;  (B)  upper  support;  (C)  control 
rod  handle;  (D)  extra  top  sheilding;  (E)  storage 
turret;  (F)  300  curie  Co-60;  (G)  shutter  shown 
open;  (H)  rear  wall;  (I)  door;  (J)  downward 
shielding;  (K)  door  carriage  (L)  door  crank;  (M) 
door  frame. 
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Fig.    3  Cross   section  of       Co  gamma  ray  source 
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Fig.    4  Decay  scheme   for  Co 
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Fig.    6  Sample  holder 
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(Hoke  catalog  No.  4611N4M).  A  rigid  aluminum  sample  holder 
(Fig.  6)  was  used  to  assure  reproducible  positioning  of 
radiolysis  vessels. 

D.  Dosimetry 

Ethylene  dosimetry  was  used  to  determine  the  absorbed 
dose  rate  of  sample  in  the  nickel  vessel.  It  has  been 
reported  that  the  hydrogen  yield  in  ethylene  under  gamma 
radiolysis  in  the  pressure  range  of  150  to  1000  torr  at 
room  temperature  is  independent  of  the  absorbed  dose  (64). 
The  G  value  for  hydrogen  in  ethelene  dosimetry  was  taken  as 
1.2  (64,65).  The  absorbed  dose  rate  in  ethylene  system  was 
determined  by  measuring  the  hydrogen  yield  at  room  tempera- 
ture, at  irradiation  times  between  4  and  24  hours,  and  at 
a  pressure  of  300  torr.  Following  irradiation,  the  amount 
of  hydrogen  (along  with  small  amounts  of  methane  and 
ethylene)  was  determined  using  a  Toepler-McLead  apparatus. 
Knowing  the  total  pressure,  volume,  temperature,  and  the 
quantities  of  methane  and  ethylene  (determined  by  gas- 
chromatographic  analysis),  the  hydrogen  yield  could  be 
calculated.  A  plot  of  the  hydrogen  yield  versus  irradia- 
tion time  gave  a  straight  line  as  shown  in  Fig.  7.  From 
the  slope  of  the  plot  and  the  accepted  G-value  for  hydrogen 
in  ethylene  (G  =  1.2)  the  absorbed  dose  rate  was 
determined.     The  absorbed  dose  rate  for  ethylene  was  2.69  x 
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Irradiation  time,  hours 

Fig.    7  Dosimetry:   hydrogen  yield   from  ethylene  as 
a  function  of   irradiation  time. 
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10-^^  eV/gh  on  December  18,  1979.  This  value  was  correc- 
ted for  Co-60  decay  during  subsequent  irradiations.  Energy 
absorption  in  propane  and  propane-oxygen  mixtures  was  cal- 
culated relative  to  the  ethylene  results.  By  application 
of  the  Bragg-Gray  cavity  principle  (66),  the  ratio  of  the 
energy  deposited  per  unit  mass  in  the  sample  to  that  in  the 
dosimeter  can  be  determined  by  the  ratio  of  these  mass 
stopping  powers  (see  appendix  I).  For  mixtures  of  propane 
and  oxygen  it  was  assumed  that  each  gas  absorbed  energy 
independently  in  proportion  to  the  mass  (more  precisely  the 
electron  stopping  power  fraction)  of  that  gas  present.  The 

dose  rates  obtained  were  1.94  x  10^^  eV/gh  for  propane  and 
19 

2.09  X  10       eV/gh  for  propane  oxygen  mixture. 


E.     Analytical  Equipment  and  Product  Analysis 

Gas  Chroma tograph  Instrument;  A  Tracer  model  1500 
gas  chromatograph  equipped  with  a  thermal  conductivity 
detector  and  flame  ionization  detector  connected  in  series 
was  used  for  the  qualitative  and  quantitative  analysis  of 
the  organic  products.  A  completely  new  injection  box 
constructed  in  this  laboratory  using  several  commercial 
components  was  fitted  onto  the  gas  chromatograph.  This  is 
shown  in  Fig.  8.  The  input  system  included  provision  for 
the  use  of  external  sample  loops,  connected  via  a  dual 
Swagelok   fitting    (obtained   from  Crawford   Fitting  Company, 
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Solon,  Ohio),  In  this  instrument,  the  gas  chromatographic 
column  is  enclosed  in  an  oven  whose  temperature  was 
controlled  by.  a  multifunction  programmer.  The  output  of 
the  gas-chromatographic  detector  was  fed  to  a  1  mV  poten- 
tiometric  recorder.  All  the  products  were  transferred  on 
the  vacuum  line  to  the  sample  loops  shown  in  Fig.  2  for 
injection  into  the  gas  chromatograph. 

Analysis  of  Products  Non-Condensible  at  -196  °C; 
After  irradiation  of  the  sample,  the  radiolysis  vessel  was 
attached  to  the  sample  loop  (S^^)  (Fig,  10)  which  in  turn 
was  attached  to  a  submanifold  leading  to  the  Toepler  purap- 
Mcleod  gauge  apparatus.  After  a  good  vacuum  had  been 
reached,  a  zero  pressure  reading  was  taken  with  an  Ealing 
Ca thetometer.  The  system  was  isolated  from  the  vacuum 
pumps  and  the  valve  of  the  nickel  vessel  was  opened.  The 
sample  loop  and  the  U-trap  were  immersed  in  liquid  nitro- 
gen. The  non-condensible  fraction  was  collected  and 
transferred  in  twelve  Toepler  pump  cycles  to  the  McLeod 
gauge  for  measurement.  The  pressure  reading  was  taken 
again.  The  gas  was  then  transferred  to  the  sample  loop 
(Fig,  9)  for  flame  ionization  gas  chromatographic  analysis. 
This  procedure  allowed  the  determination  of  small  quan- 
tities of  CH4  and  C2H4  that  contributed  to  the  pressure 
measurement.  Knowing  the  total  pressure  volume  and  tempera- 
ture, the  total  moles  of  gas  were  calculated.  Yields  of  CHy, 
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and  C2H4  were  determined  by  chromatographic  analysis.  From 
the  total  moles  of  gas  and  the  moles  of  CH4  and  C2H4,  the 
hydrogen  yield  was  determined  by  difference. 

Analysis  of  Condensible  Organic  Products;     After  the 
non-condensible  gases  were  transferred  to  the     sample  loop 
(Fig. 10),   the  liquid  nitrogen  Dewar  was  removed  from  the 
U-trap,  and  the  gases  were  transferred  to  the  sample  loop 
S-^.     The  sample  loop  S-j^  was  then  removed  from  the  vacuum 
line  and  hooked  to  the  gas  chromatograph  for  analysis.  In 
the  analysis  of  hydrocarbons,    all   the  products  were 
separated  on  a  2.7m  x  0.25  in  O.D.  stainless  steel  column 
packed  with  60/200  mesh  silica  gel,   with  a  nitrogen  carrier 
gas  whose  flow  rate  was  40  ml/min.     The  oven  compartment  in 
which  the    column    was    situated    was  at  room  temperature 
25°C.    Programming  was  started  10  minutes  after  injection 
of  the  sample  into  the  column.     The  oven  compartment  was 
programmed  at  a  rate  of  5  °C/min  up  to  a  temperature  of 
200  °C.    At  the  end  of  each  analysis,  the  column  was  condi- 
tioned at  200  °C  for  several  hours. 

Prior  to  each  analysis,  the  response  of  the  flame 
ionization  detector  was  determined  using  gaseous  propane  as 
standard.  The  flame  ionization  detector  responses  of  the 
radiolysis  products,  relative  to  propane,  are  given  in 
appendix  III.  Flame  ionization  detector  is  not  sensitive 
to  hydrogen.  It  was  measured  using  a  thermal  conductivity 
detector. 
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Fig. 10     Sample  loop   for  condensible  products. 
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In  the  identification  of  the  oxidation  products  of 
propane,  a  6m  long  carbowax-20M  (polyethylene  glycol) 
column  was  used  with  a  nitrogen  carrier  gas  flow  rate  of  40 
ml/min.  The  oven  compartment  in  which  the  column  is 
situated  was  heated  to  60  °C.  Programming  was  started  20 
minutes  after  injection  of  the  sample  into  the  column.  The 
oven  compartment  was  programmed  at  a  rate  of  4  °C/min  up 
to  a  temperature  of  200  °C.  At  the  end  of  each  analysis, 
the  column  was  conditioned  at  200  °C  for  several  hours. 

Prior  to  each  analysis,  the  response  of  the  flame 
ionization  detector  was  determined  using  acetone  as 
standard.  The  flame  ionization  detector  responses  of  the 
radiolysis  products,  relative  to  acetone,  are  given  in 
appendix  III. 

In  the  oxidation  of  propane,  a  Porpak  Q  (a  porous 
polymer  composed  of  ethyl vinylbenzene  cross-linked  with 
di viny Ibenzene  to  form  a  uniform  structure  of  a  distint 
pore  size;  type  Q  means  non  polar  phase  having  a  small  pore 
size)  column  was  also  tried,  but  the  two  major  products 
(acetone  and  propionaldehyde)  could  not  be  separated.  All 
the  radiolysis  products  were  identified  by  their  retention 
times  on  silica  gel  and  Carbowax-20M. 
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F.     Computer  Simulation 

A  computer  program  designed  to  deal  with  steady  state 
kinetics  was  used  to  test  the  validity  of  the  proposed 
mechanisms.  The  computer  simulation  program  was  originally 
written  in  FORTRAN  II  by  DeTar  (67)  and  was  modified  to 
FORTRAN  IV  and  used  on  an  S-100  bus  CP/M  based  Z80  micro- 
computer. The  FORTRAN  IV  version  of  the  program  has  been 
included  in  Appendix  IV. 

The  program  was  modified  to  create  a  disk  file  of  the 
predicted  concentration  versus  time  data,  as  well  as 
printing  it  on  the  console  video  display  or  line  printer. 
The  output  data  file  was  picked  up  from  the  disk  by  a  BASIC 
program  which  normalized  the  concentration  data  and  dis- 
played the  data  in  graphical  form  on  the  video  screen  using 
a  Vector  Graphics  high-resolution  bit-mapped  graphics  board 
on  an  Integral  Data  System  Model  225  dot-matrix  printer 
using  special  assembly  language  subroutines. 

The  input  parameters  to  the  DeTar  program  are  the 
mechanistic  steps,  their  rate  constants,  the  identification 
of  reactants  and  products,  the  initial  concentrations  and 
the  information  which  controls  the  computational  intervals. 
The  program  furnishes  the  concentration  of  each  chemical 
species  as  a  function  of  time.  The  algorithm  used  to 
perform  the  calculations  is  trapezoidal  integration  of  the 
set  of  simultaneous  rate  equations. 


III.     THE  GAMMA  RADIOLYSIS  OF  PROPANE 


A.     Experimental  Results 

The  nickel  vessels  were  filled  with  100  torr  of 
propane  gas  and  irradiations  were  carried  out  at  25  °C  in  a 
Co-60  gamma  source  over  an  absorbed  dose  range  of  1.55  x 
10^^-1.40x1021  ev  g"-^.  The  dose  rate  in  propane  was  found 
to  be  1.94x10^5  eV  g'^h'^.  Each  sample  contained  approx- 
imately 0.026  g  C3H3. 

Figure  11  shows  a  typical  chromatogram  of  the  non- 
condensible  products  and  Figure  12  shows  the  chromatogram 
of  the  condensible  organic  products.  In  all  12  products 
were  identified  by  their  retention  times.  Figures  13  to  18 
show  the  yields  of  major  and  minor  products  as  a  function 
of  dose.  These  yields  are  linear  with  absorbed  dose.  The 
major  products  are  hydrogen,  methane,  ethane,  iso-butane, 
n-butane,  iso-pentane,  n-pentane,  and  hexanes  while  the 
minor  products  are  heptanes,  octanes,  nonanes  and  decanes. 
Small  yields  of  ethylene  and  propylene  were  also  observed. 
The  G-values  of  various  products  along  with  the  results  of 
earlier  investigators  are  listed  in  Table  1. 

The  G-values  in  the  presence  of  10%  oxygen  are  listed 
in  Table  2,    along  with  the  values  reported  by  Yang  and 
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Table  1 

A  Comparison  of  the  G-Values  for  the 
Radiolysis  Products  from  Propane 


Type  of  Radiation 

Gamma 

Gamma 

Gamma 

Alpha 

Electrons 

Pressure,  cm  Hg 

10 

388 

60 

76 

4.0 

Temperature,  °C 

25 

20-25 

25 

25 

25 

Product 

This 
Work 

Birdwell 
& 

Yang 
& 

Lind 
& 

Bone , 
Blocker  & 

Crawford 

Manno 

Bardwell 

Futrell 

(1) 

(4) 

(5) 

(15) 

Hydrogen 

4.99 

(5.9)^ 

5.9 

(5.9)^ 



Methane 

1.30 

1.59 

1.48 

1.59 

1.33 

Ethylene 

trace 

trace 

0.24 



1.13 

Ethane 

1.95 

1.89 

2.01 

0.94 

1.95 

Propylene 

trace 

trace 

0.24 



1.95 

Butanes 

0.86^ 

1.00 

0.89 

1.12 

0.98 

Pentanes 

0.56*=' 

0.47 

0.41 

0.83 

0.21 

Hexanes 

0.89^ 

4.13 



1.37 

Heptanes 

0.082 

Octanes 

0.067 

Nonanes 

0.088 

Decanes 

0.033 

a.   Consisting  of 
G  =  0.25 

iso- 

butane,  G 

=  0.61 

and  n-butane. 

b.     Consisting  of 
G  =  0.14 

iso-pentane,  G 

=  0.42 

and  n-pentane. 

c.  2,3  dimethyl  butane  and  2  methyl  pentane 

d.  G-Values  based  on  the  assumption  that  G  (H,)  =  5.9 
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table  2 

G-Value  for  the  Radiolysis  Products  from  Propane 
with  Added  Scavenger. 


Products  This  Work  Yang  and  Manno 

10%  02  5%  NO. 


Hydrogen  1.87  2.48 

Methane  0.82  0.94 

Ethane  1.22  1.36 

Butanes  0.00  0.12 


Table  3 

Per  Cent  Free  Radical  Contribution 

Products  This  Work  Yang  and  Manno 

Hydrogen 
Methane 
Ethane 
Butanes 


63 
37 
17 

100 


58 
36 
32 
87 
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Table  4 

G-Values  of  Products  at  Different  Oxygen 
Concentrations  of  Oxygen  in  the  Gamma 
Radiolysis  of  Propane 


Products 


0% 


Amount  of  Added  Oxygen 


0.7% 


2% 


5% 


10% 


15% 


Hydrogen 


4.99        3.25        2.30       1.90        1.87  1.86 


Methane 


1.30        1.12        0.94       0.88        0.82  0.82 


Ethane 


1.95         1.32         1.28       1.27         1.22  1.23 
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Table  6 

Yields  of  Primary  Processes  in  Radiolysis  of  Pure 
Propane  From  Computer  Modeling  Calculation 


Process  G-value 


(1) 

{C3H8*) 

 >  C2H5*  +  CH3* 

1.00 

(2) 

C3H7*  + 

H* 

1.11 

(3) 

H2  + 

C3H6 

1.67 

(4) 

CH4  + 

C2H4 

0.83 

(5) 

C2H5*  + 

CH3- 

1.22 

(6) 

C3H5-  + 

H-     +  H2 

0.20 

(7) 

H2  + 

Polymer 

1.79 

Table  7 

Comparison  of  Measured  and  Computed  yields  in  Pure 

Propane  System 


Product 


G-Measured 


G-Computed 


Hydrogen 

Methane 

Ethane 

Butanes 

Pentanes 

Hexanes 


4.99 
1.30 
1.95 
0.86 
0.56 
0.89 


4.99 
1.32 
1.96 
0.88 
0.54 
0.89 


4A 


100.0 


to 
m 

0) 

o 
S 
o 
u 
u 


w 

•H 


80  .  0 


60.0  - 


AO.O 


20.0  - 


2.0       A.O         6.0         8.0     10.0       12.0  14.0 
Dose,   eV/g  x  10"20 

Fig.    13     Production  of  hydrogen   (pure,Z\;   10%  O2  , 
^  )    as   a  function  of   dose  in  propane 
radiolysis . 
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1  1  1  1  1  1  r 


Dose,   eV/g  x  10~20 

Fig.  14  Production  of  CH^  (pure,  O:  10%  0^,%) 
as  a  function  of  dose  in  propane  radio- 
lysis. 
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I  I 


Dose,   eV/g  x  lO"''" 

Fig,    15     Production  of  C^Hg    (pure,    □;    10%  oxygen, 
B   )    as   a  function  of  dose  in  propane 
radiolysis . 
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1  r 


0.0  3.0  6.0  9.0  12.0 

Dose,   eV/g  x  10~20 

Fig.    17     Production  of   CgH^^  A  ;    i-C3Hi2  O  a^d 

n-C5H-j^2^as  a  function  of  dose  in  propane 
radiolysis . 
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1  1  1  1  1  1  r 


Dose,   eV/g  x  lO'^O 


Fig.    18     Production  of  nonanes ,  ^  ;   heptanes, O; 

octanes  ,  O  and  decanes,Qas  a  function  of 
dose  in  propane  radiolysis. 
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Fig.    19     G-values   of  hydrogen,^;   ethane, Band 

methane, ©as  a  function  of  added  oxygen 
in  the   radiolysis   of  propane. 
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Manno  (4)  in  the  presence  of  5%  NO.  When  oxygen  is  added  as 
a  scavenger  to  propane,  it  was  found  that  the  yield  of 
hydrogen  was  reduced  by  63%,  methane  by  37%  and  ethane  by 
37%.  This  is  shown  in  Table  3.  All  the  higher  saturated 
hydrocarbons  such  as  butanes,  pentanes,  hexanes,  octanes, 
nonanes  and  decanes  were  eliminated.  Figure  19  shows  the 
variation  of  the  values  of  hydrogen  methane  and  ethane  as  a 
function  of  percent  added  oxygen. 

B.  Discussion 

Previous  reports  of  the  product  yield  in  gas  phase 
propane  radiolysis  show  considerable  disagreement  in  the 
results.  This  is  seen  in  Table  1,  where  the  results  of 
this  work  are  compared  to  the  work  of  other  investigators 
(1,4,5,15).  In  most  reported  investigations  the  material 
balance  is  not  particularly  good  with  a  hydrogen  to  carbon 
ratio  of  2.96  for  the  work  of  Bone,  Sieck,  and  Futrell 
(20),  2.82  for  Birdwell  and  Crawford  (5),  3.74  for  Yang  and 
Manno  (4),  3.71  for  Lind  and  Bardwell  (1)  and  3.14  for  this 
work. 

The  effect  of  scavenger  on  the  radiolysis  of  propane 
was  previously  reported  by  Yang  and  Manno;  they  used  NO  as 
a  scavenger  whereas  O2  was  employed  in  the  present  work. 
In  the  presence  of  oxygen,  the  hydrogen  yield  decreased  by 
63%  and  the  yields  of  methane  and  ethane  were  reduced  by 
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37%.  All  the  higher  saturated  hydrocarbons  were  elimi- 
nated, indicating  that  these  products  were  formed  from 
radicals.  From  the  results  obtained  it  is  possible  to 
propose  a  reaction  scheme  for  the  formation  of  major 
products  through  Cg. 

When  a  gas  is  irradiated,  both  excited  and  ionized 
molecules  are  formed  (Reaction  la, lb).  The  excited  propane 
can  decompose  to  give  molecular  hydrogen  and  propylene, 
methane  and  ethylene,  ethane  and  methylene,  propyl  radical 
and  hydrogen  atom,  and  an  ethyl  radical  and  a  methyl 
radical.  Reactions  2  to  6  are  well  established  through 
photochemical    studies  (20,34,73). 


^3^8 

^3^8 

(la) 

^3^8 

C3Hg+ 

+ 

e~ 

(lb) 

^3^8* 

H2 

+ 

C3H6 

(2) 

^3^8* 

CH4 

+ 

C2H4 

(3) 

C3H8* 

C2H6 

+ 

CH2 

(4) 

C3H8* 

C3H7- 

+ 

H* 

(5) 

C3H8* 

+ 

CH3* 

(6) 

Reactions  2,  3,  and  4  account  for  non-scavengable  yields 
of  H2,  CH4,  C2Hg. 

The  radicals  react  by  combination,  disproportionation 
or  abstraction  from  substrate  to  give  various  products.  Iso 
and  normal  products  are  formed  depending  on  whether  i-C3H7* 
or  n-C3H7*  takes  part  in  the  reaction.  Reactions  7,  8,  9, 
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11,  13,  and  14  have  been  proposed  by  Steacie  and  Dewar 
(73).  Reactions  10  and  12  have  been  proposed  by  Okabe  and 
McNesby  (34). 


Abstraction: 


Combination: 


+  C3H8 

> 

H2       +  C3H7' 

+  C3H8 

> 

CH4     +  C3H7* 

•  +  C3H8 

C2Hg  +  C3H7  * 

CH3* 

+  CH3* 

^2^6 

+  C2H5- 

^4^10 

CH3* 

+  C3H7- 

^4^10 

+  C3H7' 

^5^12 

^3^7' 

+  C3H7* 

^6^14 

Disproportionation; 


CH3- 

+ 

C3H7- 

CH4 

+ 

C3H6 

+ 

C2H5- 

^2^6 

+ 

C2H4 

^2^5* 

+ 

^3^7* 

C2H4 

+ 

C3H8 

^2^5* 

+ 

^3^7* 

^2^6 

+ 

C3H6 

C3H7- 

+ 

C3H7- 

C3H6 

+ 

C3H8 

C3H5. 

+ 

^2^5* 

^3^6 

+ 

C2H4 

C3H5- 

+ 

^3^7* 

C3H6 

+ 

C3«6 

(7) 
(8) 
(9) 


(10) 
(11) 
(12) 
(13) 
(14) 


(15) 
(16) 
(17) 
(18) 
(19) 
(20) 
(21) 
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When  oxygen  is  added  as  a  scavenger  it  is  seen  that 
the  yields  of  hydrogen,  methane  and  ethane  are  decreased, 
which  means  that  there  is  a  radical  yield  of  the  products. 

Products  containing  4  or  more  carbon  atoms  are 
eliminated  in  the  presence  of  oxygen  indicating  that  they 
are  all  formed  by  radical  reactions.  Reactions  7-21  account 
for  the  scavengable  yields  of  ^H^,  ^2^6'  ^4^10'  ^5^12' 
and  CgHj^^, 

Methylene  formed  in  reaction  4  inserts  itself  into 
propane  to  form  excited  butane  (Reaction  22),  which  under- 
goes decomposition  to  give  methyl  and  propyl  radicals,  and 
two  ethyl  radicals   (Reaction  23,  24). 

CH2  +  C3H8   >  C4H10*  (22) 

^4^10*   >  ^"3*      ^3^7*  (23) 

C4H10*   >  2C2H5-  (24) 

Products  containing  7  or  more  carbon  atoms  are  formed 
as  follows: 

R*  +  A   >  RA* 

2RA*   >  (RA)2 

R*   +  RA*   >  RAR 

where  R*  is  a  radical  and  A  is  an  alkene  C2H4  or  C3Hg. 

Along  with  the  radical/excited  molecular  reaction 
scheme  described  above,  numerous  ionic  processes  also 
occur,  following  primary  ionization  (Reaction  lb).  These 


55 

processes  contribute  substantially  to  certain  product 
yields.  A  sequence  of  elementary  reactions  has  been 
demonstrated  by  Bone  and  Futrell  which  lead  to  C3Hg'''  and 
CjH^'*'  as  stable  ions  (33).  In  the  mass  spectrum  of  propane 
three  ions  constitute  over  60%  of  the  total  ionization 
(18);  CjH^"^  (30.5%),  CjH^"^  (18.1%)  and  C2H3"*'  (12.3%) 

^3^8^  >   ^2^5^       ^^3*  ^25) 

 >  C2H4+  +  CH4  (26) 

 >   C2H3'''  +   H2  +  CH3.  (27) 

Hydrogen  is  also  produced  by  the  molecular  elimination 
from  C2H5"*",  C2H4'''  etc  (7).  The  various  fragment  ions 
(C2H3"'',  C2H4"*',  ^2^5^  etc.)  take  part  in  a  series  of  ion- 
molecule  reactions  involving  parent  C3H8,  which  ultimately 
lead  to  only  two  catonic  products,  C3H-7'*"  and  C3Hg'*'  (33); 
C3Hg'*"   and  C3H-7'*'  ultimately  undergo  neutralization  to  give 

^3^8  • 

C2H5+  +  C 
C2H4+  +  C 


•2^3^  +  C3H8 


8  >  ^2^6 

+  C3H7+ 

(28) 

8        ^  ^2^5*  ^3^7"*^ 

(29) 

 >  C2H6 

+  C3H6+ 

(30) 

 >  C2H4 

+  C3H/ 

(31) 

 >  C3H6 

+  C2H5+ 

(32) 

 >  C2H6 

+  C3H5+ 

(33) 

 >  C3H8 

+  C3H6+ 

(34) 

— >  C3H7- 

+  C3H/ 

(35) 
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Reactions  (25-35)  produced  portions  of  the  yields  of 
CH4,  and  C2Hg,  as  well  as  H*,  CH3*  and  C2H5*  radicals. 

Formation  of  Major  Products 

Hydrogen;  Hydrogen  is  formed  directly  from  the  de- 
composition of  excited  propane  (Reaction  2).  Molecular 
elimination  of  C3H8''',  C2H5''",  C2H4'^,  etc.  (7)  also  yields 
hydrogen.  These  reactions  account  for  the  non-scavengable 
yield  of  hydrogen,  the  G-value  of  which  is  1.87  in  the 
presence  of  oxygen. 

Hydrogen  radical  abstracts  a  hydrogen  atom  from  the 
substrate  to  give  molecular  hydrogen  (Reaction  7),  which 
accounts  for  the  scavengable  yield  of  hydrogen.  The  free 
radical  contribution  for  the  formation  of  hydrogen  is  63%. 

Methane;  Methane  is  formed  directly  by  the  decomposi- 
tion of  excited  propane  (Reaction  3).  It  is  also  formed 
from  C3H8'*"  (Reaction  26)  (7).  Reactions  3  and  26  account 
for  the  non-scavengable  yield  of  methane,  the  G-value  of 
which  is  0.82  in  the  presence  of  oxygen. 

Methyl  radicals  formed  in  reaction  6  abstract  a  hydro- 
gen from  substrate  to  give  methane  (Reaction  8)  (73).  A 
methyl  and  a  propyl  radical  disproportionate  to  yield 
methane  (15).  These  reactions  8  and  15  account  for  the 
scavengable  yield  of  methane.  The  free  radical  contribution 
for  the  formation  of  methane  is  37%. 
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Ethane:  Excited  propane  decomposes  to  give  ethane 
(Reaction  4).  It  is  also  formed  by  the  ion-molecule  reac- 
tion of  C2H5''',  C2H4'''  and  C2H3''"  with  C3H8  (Reaction  28,  30, 
and  33),  This  accounts  for  the  non-scavangable  yield  of 
ethane,  the  G-value  of  which  is  1.22  is  the  presence  of 
oxygen. 

The  ethyl  radical  formed  in  reaction  6  abstracts  a 
hydrogen  atom  from  the  substrate  to  give  ethane  (Reaction 
9)  (73).  Two  methyl  radicals  recombine  and  form  ethane 
(Reaction  10)  (34).  Two  ethyl  radicals  disproportionate  to 
give  ethane  and  ethylene  (Reaction  16).  An  ethyl  radical 
and  a  propyl  radical  can  disproportionate  to  give  ethane 
and  propylene  (Reaction  18).  All  these  reactions  6,  9,  10, 
16,  and  18  account  for  the  scavengable  yield  of  ethane. 
The  free  radical  contribution  for  the  formation  of  ethane 
is  37%. 

Butanes;     Butanes  are  formed  in  two  ways: 

(a)  Recombination  of  two  ethyl  radicals  (Reaction  11) , 

(b)  Combination  of  a  methyl  and  a  propyl  radical  (Reaction 
12)  . 

Pentanes;      Pentanes  are   formed  by  the  combination  of 
ethyl  and  propyl  radicals   (Reaction  13). 

Hexanes;      Hexanes  are   formed  by  the  recombination  of 
two  propyl  radicals   (Reaction  14). 
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Since  radicals  are  scavanged  by  oxygen,  and  butanes, 
pentanes  and  hexanes  are  not  observed  in  the  presence  of 
oxygen,  the  C4,  C5,  and  Cg  products  are  formed  by  free 
radical  reactions  (Reactions  11,  12,  13,  and  14). 

Computer  Simulation  of  the  Radiolysis  Mechanism 

A  major  goal  of  the  present  work  was  to  carry  out 
computer  modeling  of  the  reaction  kinetics,  in  order  to 
investigate  the  internal  consistency  of  the  yield  data,  as 
well  as  the  adequacy  of  a  simple  reaction  scheme  to 
describe  the  data.  The  mechanism  consisting  of  reactions 
given  in  Table  5  was  tested  using  rate  constants  given  in 
the  literature.  Since  the  entire  scheme  is  too  complex  to 
be  suitable  for  modeling,  several  simplifications  were  made 
in  selecting  a  subset  for  the  computer  work.  In  parti- 
cular, it  was  assumed  that  ion-molecule  processes  were 
completed  before  onset  of  the  secondary  radical  reaction 
stage,  and  that  the  ionic  processes  in  effect  constitute 
additional  pathways  to  the  same  set  of  primary  products  and 
radical  intermediates  as  formed  via  neutral  fragmentation 
and  elimination  reactions.  (It  was  necessary  to  account 
explicitly  for  the  role  of  C2H5+  in  the  mechanism,  as 
explained  below.) 
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Formation  of  and  C2  radicals  is  represented  as 
occurring  via  a  two  step  sequence  (Reaction  SI  and  S7  in 
Table  5)  due  to  the  observation  that  essentially  no  oxygen 
containing  C-^  and  C2  products  are  observed  in  the  presence 
of  added  oxygen.  The  species  designated  C2H5*  in  step  5  is 
in  fact  known  to  be  C2H5"*'  ion,  which  is  a  precursor  to  a 
(nearly)  unsca vengable  C2Hg  yield  via  an  H"  transfer  pro- 
cess. Reaction  S8  is  a  simplified  representation  of  this 
process;  it  ignores  charge  state  and  the  problem  of  ion- 
neutralization.  It  is  thus  assumed  that  the  C^^-j"*"  ion 
actually  formed  in  Reaction  S8  is  a  precursor  to  C3H-7* 
radical  via  an  unspecified  reaction  path.  Reaction  S6  has 
been  added  to  the  scheme,  because  allyl  radical  plays  an 
important  role  in  the  radiolytic  oxidation  of  propane,  in 
the  formation  of  allyl  alcohol  and  acraldehyde. 

Reactions  S1-S6  were  assigned  pseudo-rate  constants  on 
the  bases  of  the  G  values  given  in  Table  5.  That  is,  the  G 
value  was  converted  to  an  actual  rate  of  product  formation 
in  units  of  moles,  liters,  and  seconds,  using  the  known 
dose  rate  in  the  system. 

The  apparent  first  order  rate  constants  (kpsg^jj^)  for 
the  primary  reactions  are  calculated  as  follows: 

(Molecules/sec)  x  100 

G  =  

(eV  deposited/liter  sec) 
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Molecules/sec  = 


'^Pseudo 


G  X  (ev/liter  sec) 


100 


Moles         G  X   (eV/liter  sec) 


liter  sec         100  x 


Moles 


Let  Y  =    =  Rate  =  K  Pseudo  [C3H8] 

liter  sec 


Rate 


[C3H8]  [C3H8] 


If  the  approximate  G-value  is  known,  Y  which 
is  in  moles/(liter  sec)  can  be  calculated;  and  knowing  the 
concentration  of  C3H8  in  moles/lit,  kpg^^^^^  can  be  calcu- 
lated. 

This  procedure  was  required  since  the  kinetic  integra- 
tor package  can  handle  only  first  or  second  order 
processes.  Since  the  propane  concentration  is  virtually 
constant,  the  procedure  described  leads  to  a  constant  rate 
of  formation  of  the  species  in  question.  (Third  order 
processes  usually  involve  the  substrate  as  a  third  body; 
they  are  conveniently  treated  as  second  order  by 
multiplying  the  third  body  concentration  into  the 
termolecular  rate  constant,  giving  an  effective  bimolecular 
rate  constant.) 


61 


The  rate  constant  for  Reaction  S7  was  given  the  value 
of  9,0x10^  sec""^  in  the  present  study,  since  further  work 
on  oxygen-containing  products  suggested  that  this  reaction 
may  be  subject  to  interference  (via  collisional  deactiva- 
tion of  C3Hg*)  at  high  concentrations  of  added  This 
process  will  be  the  exclusive  fate  of  C3Hg*  in  the  unscav- 
enged  system,  since  no  competitive  fate  of  the  reactant  is 
indicated.  The  radical  products  would  presumably  be 
trapped  with  any  substantial  concentration  of  added  oxygen. 
Reaction  S8  was  assigned  an  appropriate  high  rate  constant 
for  an  efficient  ion-molecule  reaction,  and  Reactions  S9 
and  S10  were  given  arbitrary  rate  constants  of  1x10^  liter 
mole"^  sec~^.  Rate  constants  in  this  range  permit  inter- 
ference by  oxygen  in  the  1  to  15%  concentration  range,  as 
observed  in  the  radiolytic  oxidation  of  propane. 

Reactions  Sll  to  S21  were,  in  general,  assigned  the 
best  available  values  obtained  from  the  literature  (68-72). 
There  was  a  problem  in  the  case  of  reaction  S16,  however, 
since  published  rate  constants  varied  over  a  wide  range 
from  2.24x10^  liter  mole"^  sec"^  to  6.6x10^^  liter  mole"-^ 
sec~^;  we  used  an  intermediate  values  of  1.4x10^0.  Reaction 
S19  has  published  values  which  vary  from  3.16x10^  to 
6.31x10^  liter  mole"^  sec"^.  We  tried  several  values  which 
lie  in  the  range  of  1.59x109  liter  mole"^  sec~^  to  6.31x10^ 
liter  mole"!  sec"!  (72).  the  value  used  in  computing  the 
curves  in  Figure  20  and  21  was  3.4x10^  liter  mole"!  sec"!. 
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We  also  assigned  arbitrary  values  of  1x10"^  to  Reactions  S12 
and  S13  because  the  literature  values  (28  and  1.48  liter 
mole~l  sec"^)  appeared  implausibly  small.  In  fact,  use  of 
these  literature  values  gave  computed  yield  curves  in  clear 
disagreement  with  the  experimental  results  of  the  present 
work . 

The  input  to  the  computer  program  was,  essentially, 
the  yields   of   molecular   hydrocarbon  products    from  the 
oxygen  scavenged  system  plus  the  set  of  thermal  radical 
Reactions  S12-S21,   along  with  their  rate  constants.  Yields 
for  two  primary  processes  leading  to  radical  intermediates 
(Reaction  SI  and  S2)  were  varied  to  achieve  the  best  fit  of 
product  yields    in   the   absence  of  scavenger.    About  100 
computer  runs  were  carried  out.  The  yields  for  the  primary 
processes  used  in  the  best-fit  computer  calculations  are 
given  in  Table  6.  Table  7     shows  a  comparison  of  measured 
and  computed  yields  of  the  products.     Figures  20  and  21 
show  the  plot  of  concentration  versus  time  of  the  various 
products  obtained  by  computer  simulation.    A  print  out  of 
the  mechanism,      their  rate  constants,    concentration  of 
products  at  various  time  periods  and  the  data  sheet  are 
given  in    Appendix  V. 

Although  for  processes  S7-S10  rate  constants  are 
listed  as  "arbitrary"  in  Table  5,  these  parameters  are  not 
"degrees  of  freedom,"  subject  to  adjustment  for  improvement 
of  the  predicted  results.  The  reason  is  that  Reactions  S7, 
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Fig.  20 


Computer  predicted  yields  of  hydrogen, A; 
ethane,  □  and  hexanes  ,  ^  as  a  function  of 
radiolysis  time. 
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Fig.    21     Computer  predicted  yields   of  methane,  O; 

butanes,  0  and  pentanes, <^  as  a  function  of 
radiolysis  time. 
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S8,  and  the  pair  S9-S10  are  the  only  reaction  channels 
specified  for  C3Hg*,  C2H5*,  and  H*  respectively.  The  only 
adjustment  which  would  affect  predictions  is  the  ratio 
kS9/kS10.  Since  yields  of  five  products  (excluding 
hydrogen)  were  predicted  with  only  two  entirely  free  vari- 
ables (Gg2  and  Gq2^ '  clear  that  the  computer  integra- 
tion computation  not  only  tests  consistency  between  the 
reaction  scheme  and  the  predicted  yields,  but  also  has  at 
least  some  predictive  value,  (There  was  some  adjustment 
made  to  other  rate  conatants,   as  mentioned  previously.) 

The  situation  regarding  the  hydrogen  yield  requires 
separate  discussion.  In  this  study  and  essentially  in  all 
other  investigations  of  the  gas  phase  radiolysis  of 
saturated  hydrocarbons,  there  is  a  material  balance 
problem,  in  that  the  hydrogen  yield  is  considerably  in 
excess  of  the  combined  yields  of  unsaturated  and  telomeric 
hydrocarbons  needed  to  account  for  it  on  the  basis  of 
stoichiometry . 

This  problem  has  been  disposed  of  arbitrarily  in  the 
modeling  calculations  by  postulating  the  deposition  of 
olefins  on  the  vessel  walls  (presumably  as  polymer). 
Although  arbitrary,  this  assumption  may  well  be  essentially 
correct,  since  strong  heating  of  the  radiolysis  vessels 
after  removal  of  volatile  products  usually  produces  consi- 
derable outgassing,  presumably  of  pyrolysis  products. 
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An  essentially  similar  postulate  was  presented  earlier 
by  Bone,  Sieck,  and  Futrell  (19),  who  concluded  that  a  sub- 
stantial part  of  the  hydrogen  yield  arises  from  neutraliza- 
tion of  the  "unreactive"  ion  C^^-j"*^  at  the  vessel  walls. 
They  did  not  discuss  in  detail  what  becomes  of  the  residual 
hydrocarbon  moieties,  but  it  appears  that  most  of  it  ends 
up  as  polymer  on  the  walls,  on  the  basis  of  the  material 
balance  situation. 

As  mentioned  above,  the  ionic  processes  implied  by  but 
not  explicitly  included  in  our  reaction  scheme  have  been 
discussed  in  detail  elsewhere  (19,29-33).  It  is  also 
known,  particularly  from  the  work  of  Ausloos  and  coworkers, 
that  there  is  an  additional  neutral  fragmentation  process 
which  we  did  not  include 

^3^8     — — — ^  ^2^6  ^^2 

The  ethane  formed  is  a  part  of  the  "molecular"  yield 
and  the  carbene  inserts  into  substrate,  giving  a  "hot" 
butane  which  can  be  stabilized  or  fragmented  into  two  ethyl 
radicals  or  a  methyl  plus  a  propyl  radical.  Some  computa- 
tions were  performed  including  these  steps,  but  ultimately 
it  was  decided  to  omit  them  for  the  practical  reason  that 
the  complexity  of  the  scheme  was  beginning  to  tax  the 
capabilities  of  the  relatively  simple  kinetic  integrator. 
For  the  same  reason,  no  distinction  was  made  between  n- 
propyl   and   i-propyl   radicals   in  the   radiolysis  of  pure 
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propane,  (The  latter  distinction  came  to  be  of  major  impor- 
tance in  the  oxygen-containing  products  in  the  radiolytic 
oxidation  of  propane,  since  i-propyl  radical  is  precursor 
to  i-propyl  alcohol  and  acetone,  whereas  n-propyl  radical 
forms  n-propyl  alcohol  and  propionaldehyde.) 

Effectively,  the  kinetic  integrations  performed  in  the 
computation  scheme  concern  reactions  occurring  on  the  10"^ 
to  10  s  time  scale,  integrated  over  minutes  to  hours.  In 
principle,  it  would  be  feasible  to  carry  out  computations 
on  the  10~12        10"^  s  time  scale,  which  would  use  mass 

spectrometric  fragmentation  patterns  and  W  values  as  input 
data  (along  with  information  from  VUV  photochemistry  for 
the  neutral  processes)  in  conjunction  with  ion-molecule 
reaction  rate  data.  An  attempt  to  do  this  was  not  made, 
mainly  because  most  of  the  ion-molecule  reaction  work  on 
propane  was  done  in  the  rather  qualitative  1960  -  1965  era; 
for  the  most  part,  reaction  pathways  were  identified,  but 
rate  constants  (or  cross  sections)  were  not  measured. 

C.  Summary 

The  gamma-radiolysis  of  gaseous  propane  was  studied  at 
100  torr  pressure  and  25  °C,  both  pure  and  with  10%  added 
oxygen.  For  the  pure  system,  the  major  radiolytic  products 
and  their  respective  G-values  were  hydrogen,   4.99;  methane. 
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1.30;  ethane,  1.95;  i-butane,  0.61;  n-butane,  0.25;  i- 
pentane,  0.42;  n-pentane,  0.14  and  hexane,  0.89.  The  minor 
products  included  heptanes,  octanes,  nonanes  and  decanes. 
Small  yields  of  ethylene  and  propylene  were  also  observed. 
Several  of  the  higher  products  had  not  been  reported  pre- 
viously. It  was  observed  that  the  addition  of  10%  oxygen 
reduced  the  hydrogen  yield  by  63%,  methane  by  37%  and 
ethane  by  37%.  All  higher  saturated  hydrocarbons  were 
eliminated  in  the  presence  of  oxygen.  The  results  were 
discussed  in  terms  of  both  excited  and  ionized  propane.  In 
the  mass  spectrum  of  propane  three  ions  constitute  over  60% 
of  the  total  ionization,  C2H5''' (30.5%)  ,  02^/^'^  (18.1%)  and 
C2H3"*'  (12.3%).  The  fragment  ions  then  undergo  ion  molecule 
reaction  with  propane  which  contribute  substantially  to 
certain  product  yields.  The  non-scavengable  products  were 
accounted  for  primarily  by  the  direct  interaction  of 
propane  with  gamma  rays  to  give  molecular  products  and  by 
various  ion-molecule  reactions.  The  scavangable  products 
were  attributed  to  the  radical  fragments,  which  react  by 
combination,  disproportionation  or  abstraction  from  the 
substrate.  The  reaction  scheme  for  the  formation  of  major 
products  through  Cg  was  examined  using  computer  modeling 
based  on  a  mechanism  involving  24  reactions  and  their  rate 
constants.  Yields  could  be  brought  into  agreement  with  the 
data  within  experimental  error  in  nearly  all  the  cases,  but 
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it  was  necessary  to  assume  that  the  molecular  hydrogen 
yield  was  accompanied  by  the  deposition  of  polymer  on  the 
vessel  walls. 


IV.     RADIOLYTIC  OXIDATION  OF  PROPANE 


A.     Experimental  Results 

At  the  onset  of  research  there  was  no  preconceived 
idea  what  type  of  oxygenated  products  would  be  formed  in 
the  radiolytic  oxidation  of  propane;  therefore  a  wide  range 
of  standards,  which  included  oxygenated  hydrocarbons  such 
as  ketones,  aldehydes,  alcohols,  acids,  esters,  ethers, 
etc.  were  injected  into  the  gas-chromatograph.  Carbowax- 
20M  and  Porapak  Q  columns  were  used  with  a  nitrogen  carrier 
gas  flow  rate  of  40  ml/min.  The  temperature  of  the  oven 
was  initially  set  at  60  °C  for  Carbowax-20M  and  at  80  °C 
for  Porapak  Q.  Programming  was  started  10  minutes  or  20 
minutes  after  injection,  for  Porapak  Q  and  Carbowax-20M, 
respectively.  This  procedure  gave  the  best  separation  of 
compounds.  Retention  times  and  retention  temperatures  for 
various  standards  were  noted  as  given  in  Table  8. 

After  the  propane-oxygen  mixtures  were  irradiated  for 
different  time  periods,  the  products  were  transferred  to  a 
sample  loop  on  the  vacuum  line;  the  sample  loop  was  then 
attached  to  the  gas-chromatograph  for  analysis.  The  pro- 
ducts were  identified  by  comparing  their  retention  time  and 
retention  temperatures  with  the  standards.     The  products 
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identified  in  the  radiolytic  oxidation  of  propane  (in- 
dicated by  asterisks  in  Table  8)  were  mainly  ketones, 
aldehydes,  and  alcohols,  which  were  similar  to  the  products 
formed  in  the  autoxidation  of  hydrocarbons.  Intially  both 
the  above  mentioned  columns  were  used  in  the  identification 
of  products,  but  only  the  Carbowax-20M  column  was  used  for 
the  quantitative  determination  of  products  because  acetone 
and  propionaldehyde  (which  were  major  products)  were 
separable  on  Carbowax-20M  but  not  on  Porapak  Q.  Figure  22  a 
and  b  show  a  chromatogram  of  the  condensible  organic 
products. 

The  nickel  vessels  were  filled  with  a  mixture  of 
propane  (100  torr)  and  oxygen  (10  torr)  and  irradiations 
were  carried  out  at  26  °C  in  a  Co-60  gamma  source  over  an 
absorbed  dose  range  of  2.5x1020-  15.22x1020  eV/g.  The 
dose  rate  was  found  to  be  2.11x1019  eV/gh.  Each  sample 
contained  approximately  0.028  g  of  mixture  (propane  and 
oxygen) , 

In  all  11  oxygen  containing  products  were  idenfied  by 
their  retention  time  and  temperatures.  Figures  23-27  show 
the  yields  of  major  and  minor  products  as  a  function  of 
dose,  these  yields  are  linear  with  absorbed  dose.  The 
major  products  are  acetone,  iso-propyl  alcohol,  propion- 
aldehyde, acraldehyde,  n-propyl  alcohol  and  allyl  alcohol, 
while  the  minor  products  are  iso-butyl  alcohol,  tert-amyl 
alcohol,   n-butyl   alcohol,    iso-amyl   alcohol,   and  n-amyl 
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Table  8 

Retention  Times  and  Retention  Temperatures 
for  Various  Standards 


Column 

Carbowax-20M 

Porapak  Q 

Compound 

Retention 

1  X  III  c  ^  111  X  11 

Retention 

o  m             0  ^ 
i  c  ni  p  •  /  u 

Retention 

XxIUC/  luiri 

Retention 

'P  ^  m  O 
1  6111  p*  , 

Formaldehyde 

3.25 

60 

9.14 

80 

Acetaldehyde 

4.30 

60 

12.51 

115 

Propionaldehyde* 

7.10 

60 

20.40 

132 

I-Butyraldehyde 

10  .  50 

60 

26 .10 

158 

TO     1  C 

27  #  45 

lo2 

n— Valeraldehvde 

•51    yi  c 

j2  •  20 

193 

Acraldehvde* 

1  a  ODi 

DM 

128 

Acetone* 

o  .  Z3 

C  fli 

20  .  40 

132 

Methyl-Ethyl 
ketone 

12.50 

60 

27.45 

162 

1 J  .  ^19 

O0 

9.35 

80 

E thanol 

1  C     T  fli 

C  ft 

60 

16.50 

115 

I— ProDvl  Alcohol* 

14.30 

60 

21-45 

N-Propyl  Alcohol* 

26.45 

90 

24.15 

149 

Sec-Butyl 
Alcohol 

26.00 

89 

23.30 

146 

I-Butyl  Alcohol* 

31.05 

110 

29.25 

172 

N-Butyl  Alcohol* 

34.30 

128 

30.40 

180 

Tert-Amyl* 
Alcohol 

24.45 

80 

34.00 

198 
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Table  8  (continued) 


Column 

Carbowax-20M 

Porapak  Q 

Compound 

Time,min 

Temp.,  °C 

Ke uenuion 
Time,  min 

Ke tent 1 on 
Temp.,  °C 

I-Amyl  Alcohol* 

Q  701 

X'kO 

io  »  iVi 

O  U  1 

202 

N-Amyl  Alcohol* 

AOi  1  a 

HV  »  ±v 

^  Ati 

1  fl  o 

ob  .  3D 

210 

I-Hexyl  Alcohol* 

X  O  3 

JO  .  ZD 

210 

N-Hexyl  Alcohol* 

n't » Vt) 

1/10 

A  ^  fxn 

210 

Allyl  Alcohol* 

■J  ^  m  HV 

X  i.  o 

O  /I  CI 

Z4  .31 

151 

Propylene  oxide 

a  a 
K>V) 

1  o  cm 
l7  .  30 

130 

Formic  Acid 

"^Q  Aa 

1  ^  Q 

1  /  .  DD 

120 

Acetic  Acid 

aa 

1  n  CI 

1 0  19 

Z  J  .30 

1  il  C 

145 

Propionic  Acid 

AO  aa 

l7  J 

30.00 

176 

Methyl  Acetate 

7  •  X  3 

23 . 30 

145 

Ethyl  Acetate 

J. Z  .  3D 

28  .  30 

167 

I-Propyl  Acetate 

14.00 

60 

32.10 

190 

N-Propyl  Acetate 

21.50 

70 

34.00 

198 

N-Butyl  Acetate 

30.45 

110 

Ethyl  Formate 

9.30 

60 

N-propyl  Formate 

15.45 

60 

Methyl 
Propionate 

13.55 

60 

Ethyl 

Propionate 

19.45 

60 
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Table  8  (continued) 


Column 

Carbowax-20M                  Porapak  Q 

Compound 

Retention 
Time,min 

Retention    Retention  Retention 
Temp.,  °C  Time,   min     Temp.,  °C 

N-propyl 
Propionate 

28.45 

100 

Di-Ethyl  Ether 

3.45 

60 

I-Propyl  Ether 

4.20 

60 

N-Propyl  Ether 

4.25 

60 

N-Butyl  Ether 

19.30 

60 

Allyl  Ether 

4.25 

60 

*  Identified  in  the  radiolytic  oxidation  of  propane. 
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Table  9 

G-values  for  the  Products  Containing  Oxygen 
in  the  Radiolytic  Oxidation  of  Propane 
When  10%  Oxygen  Is  Added 


Products  G~values 


Acetone 

0.98 

Iso-Propyl  Alcohol 

0.86 

Propionaldehyde 

0.43 

N-Propyl  Alcohol 

0.11 

Acralaldehyde 

0.14 

Allyl  Alcohol 

0.038 

Iso-Butyl  Alcohol 

8.5  X 

10" 

■3 

Tert-Amyl  Alcohol 

7.4  X 

10" 

•3 

N-Butyl  Alcohol 

2.9  X 

10" 

■3 

N-Amyl  Alcohol 

1.6  X 

10" 

3 
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Table  11 

Experimental  G-Values  of  Products  Containing  Oxygen 
at  Different  Oxygen  Concentration 


Amount  of  Added  Oxygen 

Products  _—  

1%  3%  5%  7%  10%        13%  15% 


Acetone        1.20      1.11      1,07      1.04      0.98      0.97  0.92 


I-Propyl       1.11       1.06       0.96       0.93       0.86       0.84  0.80 

Alcohol 


Propion-      0.33       0.38       0.40       0.42      0.43       0.43  0.44 

Aldehyde 


N-Propyl      0.16      0.14      0.13      0.12      0.11      0.11  0.11 

Alcohol 


Acral-  0.13      0.14      0.14      0.15      0.15      0.15  0.15 

dehyde 


Allyl  0.044     0.044     0.041     0.038     0.038     0.037  0.036 

Alcohol 
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Table  12 

Computer  Predicted  G-Values  of  the  Products  Containing 
Oxygen  at  Different  Oxygen  Concentration 


Amount  of  Added  Oxygen 
Products  

1*  3%  5%  7%  10%  13%  15% 

Acetone      1.13      1.09      1.05      1.02        0.98        0.96  0.94 


I-Propyl     1.06       1.01       0.97       0.94        0.90  0.88 

Alcohol 


0.86 


Propion-     0.31       0.35       0.37       0.40        0.43        0.46  0.47 

Aldehyde 


N-Propyl     0.14      0.13      0.13      0.12        0.12        0.11  0.11 

Alcohol 


dehjdi        ^'^^      ^'^^      ^'^^      ^'^^        ^'^^        ^-^^  ^-^5 


Allyl  0.042   0.041   0.039     0.038   0.037     0.036  0.036 

Alcohol 


Table  13 

Yield  of  Primary  Processes  in  the  Radiolytic  Oxidation  of 
Propane  From  Computer  Modeling  Calculation 


Process  G-value 


(1)  C3H8  >  C3H8*  1.00 

(2)  C3H8  >  i-C3H7'  +  H*  1.55 

(3)  C3H8  >  n-C3H7*  +  H*  0.47 

(4)  C3H8  >  C3H5'       +  H»  +  H2  0.20 


Table  14 

Comparison  of  Measured  and  Computed  Yields  in  the 
Radiolytic  Oxidation  of  Propane 


Product  G-Measured  G-Calculated 

Acetone  0.98  0.98 

I-Propyl  Alcohol  0.86  0.90 

N-Propyl  Alcohol  0.11  0.12 

Propionaldehyde  0.43  0.43 

Acraldehyde  0.15  0.15 

Allyl  Alcohol  0.038  0.037 
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Fig.    23     Production   of  acetone,  O  and  i-propyl 

alcohol, I  as  a  function  of  dose  in  the 
radiolytic  oxidation  of  propane. 
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Fig.  24 


Production  of  propionaldehyde  D  and 
n-propyl  alcohol^  as   a  function  of  dose 
in  the  radiolytic  oxidation  of  propane. 
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T  1  1  1  1  r 


2.0       4.0       6.0       8.0     10.0     12,0  14.0 

Dose,   eV/g  x  10~^° 

Fig.    25     Production  of  acraldehyde  A  and  allyl 
alcohol  <0>  as   a   function  of  dose  in  the 
radiolytic  oxidation  of  propane. 
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Fig, 
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26     Production  of  i-butyl  alcohol,  O  and  tert- 
amyl  alcohol,  □  as   a  function  of  dose  in 
the  radiolytic  oxidation  of  propane. 
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T  1  1  1  r 


2.0       4.0       6.0       8.0       10.0     12.0  14.0 
Dose,   eV/g  x  10~20 

Fig.    27     Production  of  n-butyl  alcohol,  B;  n-amyl 
alcohol ,®  and  i-amyl  alcohol^as  a 
function  of  dose  in   the   radiolytic  oxida- 
tion of  propane. 
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alcohol.  Small  yields  of  iso-hexyl  alcohol  and  n-  hexyl 
alcohol  were  also  observed.  The  G-values  of  various 
products  at  10%  added  oxygen  are  listed  in  Table  9.  There 
was  no  apparent  change  in  the  G-values  at  pressures  of  50, 
100,  150  torr.  When  the  concentration  of  oxygen  is  varied, 
it  is  seen  that  at  lower  concentrations  of  oxygen  the  yield 
of  acetone,  iso-propyl  alcohol,  and  n-propyl  alcohol 
increases,  the  yield  of  propionaldehyde  decreases  while  the 
yield  of  acraldehyde  and  allyl  alcohol  remains  the  same. 
This  is  shown  in  Figs.  30a  and  31a.  The  G-values  of 
various  major  oxygen  containing  products  at  different  con- 
centrations of  oxygen  at  constant  dose  10.1x10^^  eV/g  are 
listed  in  Table  11. 

B.  Discussion 

If  one  takes  into  account  all  ionic  and  free  radical 
reaction  pathways  known  to  take  part  in  the  gamma  radio- 
lysis  of  gas  phase  propane,  the  system  would  appear  to  be  a 
rather  complex  one.  It  might  seem  that  addition  of  to 
the  system  would  lead  to  even  greater  complexities.  In  the 
early  stages  of  interpretation  of  the  data,  however,  it  was 
discovered  that  the  product  spectrum  and  even  relative 
product  yields  bear  a  remarkable  resemblance  to  the  species 
formed  in  the  low  temperatures  autoxidation  of  gaseous 
propane.    Accordingly,  an  attempt  was  made  to  interpret  the 
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data  on  the  basis  of  a  simple  free  radical  mechanism.  As 
will  be  seen  below  this  effort  was  quite  successful  in 
providing  both  a  qualitative  and  quantitative  description 
of  the  product  spectrum  and  experimental  G-values  in  the 
radiolysis  of  propane  oxygen  mixtures. 

It  is  reasonable  to  inquire  what  is  the  physical  basis 
for  the  apparent  simplicity  of  this  system;  i.e.  why  is  it 
not  necessary  to  consider  ionic  reaction  pathways  in  inter- 
preting the  results?  The  most  obvious  interpretation  is  to 
suggest  that  O2  acts  in  the  system  not  only  to  intercept 
thermal  free  radicals,  but  also  to  capture  thermal  elec- 
trons. In  consequence,  the  high  energy  neutralization  path- 
ways which  must  be  present  in  the  pure  system 

^3^8^  +  e~   >    products  (a) 

are   replaced  with  a  new  neutralization  pathway 

^3^8^  +   >  products  (b) 

The  latter  process  is  characterized  by  an  energy 
release  which  is  lower  to  the  extent  of  the  electron 
capture  energy  of  Oj.  Several  consequences  would  follow; 
any  slow  fragmenting  processes  involving  long-lived  excited 
states  of  CjHg"^  would  be  quenched,  and  fragmentation  pro- 
cessess  resulting  from  the  relatively  high  energy  neutrali- 
zation process  (a)  would  be  replaced  by  new  reaction  path- 
ways subsiquent  to  the  lower  energy  neutralization  process 
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(b).  Not  only  is  the  total  energy  release  in  process  (b) 
less  than  process  (a),  but  the  exothermici ty  might  be 
shared  by  the  C^Hq'*'  and  O2"  species  involved,  and  in  fact 
the  reaction  might  involve  H"*"  transfer  to  Oj"/  with  the 
consequence  that  essentially  all  the  energy  released 
involves  (or  neutral  H02*). 

It  seems  clear  that  addition  of  serves  to  quench 
many  ionic  processes  in  the  system.  This  includes  not  only 
the  prevention  of  subsequent  fragmentation  of  longer-lived 
states  of  C3Hg"*',  but  also  the  prevention  of  subsequent  ion- 
molecule  Reaction  on  the  part  of  first  generation  fragmen- 
tation ions  formed  from  C3Hg"''.  While  sufficient  energy 
might  be  released  in  neutralization  step  (b)  to  cause 
fragmentation  of  CjHg*,  the  available  energy  may  be  suffi- 
ciently low  to  prevent  extensive  secondary  fragmentation. 
In  fact,  neutralization  via  the  H"*"  transfer  route  appears 
especially  attractive,  since  this  would  yield  exclusively 
species  with  the  stoichiometry  C3H-7*. 

There  is  another  favorable  aspect  to  the  model  just 
presented.  As  described  below  in  connection  with  computer 
simulation  work,  it  was  observed  that  the  total  G-value  for 
all  processes  modifying  C3H8  is  approximately  the  same  in 
the  scavenged  and  unscavenged  system.  However,  it  was  found 
that  the  detailed  primary  yields  from  the  unscavenged 
system  could  not  be  carried  over  without  modification  to 
the  oxygen  scavenged  system.  An  additional  aspect  of  the 
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situation  is  that  among  the  processes  wiped  out  by  adding 
oxygen  appear  to  be  those  processes  which  in  the 
unscavenged  system  lead  to  C-C  bond  rupture.  This  observa- 
tion is  perhaps  the  best  evidence  that  the  actual  mechanism 
of  neutralization  step  (b)  is  the  H"*"  transfer  since  as 
noted  above  this  would  give  only  C2ii'j'  radical,  and  would 
not  necessarily  give  a  large  deposition  of  energy  in  the 
hydrocarbon.  In  fact,  the  neutralization  would  occur  in  the 
(H"^  ...  02~)  reaction  complex. 

Formation  of  Products  Via  Free  Radical  Mechanism 

Oxygen  readily  adds  to  alkyl  free  radicals,  to  give 
relatively  stable  peroxy  radicals  (36,61). 

R*  +  O2   >  R-0-0' 

From  the  results  in  Table  9  it  is  clear  that  three 
carbon  species  predominate  among  the  products;  all  the 
species  formed  with  G-values  greater  than  0.01  contain 
three  carbon  atoms.  These  species  formed  can  be  derived 
from  i-propyl,   n-propyl,   and  allyl  radicals. 

C3H8  >  C3H8*  (1) 

C3H8  — — ->  i-C^H^'  +  H*  (2) 

C3H8  >  n-C3H7-  +  H*  (3) 

C3H8  >  C3H5-       +  H2  +  H-  (4) 
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The  formula  C3H8*  represents  a  rather  long-lived 
excited  state  of  C3H8/  which  is  subject  to  Reaction  with  O2 
(or  collisional  quenching  by  O2);  the  role  of  this  species 
is  discussed  below. 

Reactions  2  and  3  are  shown  as  producing  hot  H*  ex- 
clusively, along  with  n-  and  i-propyl  radicals.  It  is  not 
critical  that  some  fraction  of  the  hydrogen  produced  may  be 
thermal,  since  a  portion  of  the  H*  is  thermalized  in  any 
event.  (See  Reaction  21  and  related  discussion  given 
below).  Reaction  4  seems  to  be  the  most  straightforward 
process  leading  to  the  formation  of  allyl  radical,  C3H5*, 
which  is  precursor  to  two  observed  products  of  moderate 
importance. 

These  radicals  i-C3H7*,  n-C3H7*  and  C3H5*  formed  in 
Reactions  2,  3,  and  4  combine  with  oxygen  to  give  i-propyl 
peroxy,   n-propyl  peroxy  and  allyl  peroxy  radicals. 


i-C3H7- 

+ 

O2 

i-C3H7-0-0* 

(5) 

n-C3H7* 

+ 

O2 

n-C3H7-0-0* 

(6) 

C3H5- 

+ 

O2 

C3H5-0-0* 

(7) 

Alkyl  peroxy  radicals  combine  pairwise  to  give  alkoxy 
radicals  along  with  molecular  oxygen  (36). 

2  R02*   >  2  RO'  +  O2 
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For  example,  two  i-propyl  peroxy  radicals  give  iso- 
propoxy  radicals  plus  molecular  oxygen,  n-propyl  peroxy 
gives  n-propoxy  and  allyl  peroxy  gives  allyloxy  radical 
products. 

2  i-C3H7-0-0*   >  2  i-C3H7-0'  +  '  (8) 

2  n-C3H7-0-0*   >  2  n-C3H7-0'  +  O2  (9) 

2  C3H5-0-0-   >  2  C3H5-0-       +  (10) 

The  cross-reactions  must  also  occur:  i-C3H7-0-0*  can 
combine  with  n-C3H7-0-0-  or  C3H5-0-0-;  and  n-C3H7-0-0*  can 
combine  with  C3H5-0-0*  to  give  molecular  oxygen  and  the 
corresponding  alkoxy  radicals. 

Combination  of  two  alkoxy  radicals  to  give  a  peroxide 
appears  less  probable  than  disproportionation  to  give  an 
alcohol  and  a  carbonyl  compound  (36). 

2  RCHjO*   >  R-CH2-O-O-CH2-R  less  probable 

2  RCH20-   >  RCH2OH  +  RCHO  more  probable 

When  two  i-propoxy  radicals  combine  the  products  are 
acetone  and  iso-propyl  alcohol. 

2   i-C3H70-  >    (CH3)2C=0       +         i-C3H70H  (H) 

acetone  i-propyl  alcohol 

Two  n-propoxy  radicals  combine  to  give  propionaldehyde 
and  n-propyl  alcohol. 
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2  n-C3H70*   >  CH3CH2CHO      +        n-C3H70H  (12) 

propion-         n-propyl  alcohol 
aldehyde 

When  an  i-propoxy  radical  combines  with  a  n-propoxy 
radical  the  products  are  either  (a)  acetone  and  n-propyl 
alcohol,   or  (b)  propionaldehyde  and  iso-propyl  alcohol. 

i-C3H70*  +  n-C3H70-   >   (CH3)2C=0      +      n-C3H70H  (13) 

acetone  n-propyl  alcohol 

i-C3H70*  +  n-C3H70-   >   i-C3H70H        +       CH3CH2CHO  (14) 

iso-propyl  Propion- 
alcohol  aldehyde 

An  i-propoxy  radical  can  combine  with  an  allyloxy 
radical  in  two  ways  to  give 

a)  acetone  and  allyl  alcohol  or  b)  i-propyl  alcohol  and 
acraldehyde. 

i-C3H70-  +  C3H50-   >    (CH3)2C=0     +     CH2=CHCH20H  (15) 

acetone  allyl  alcohol 

i-C3H70-  +  C3H50-   >  I-C3H7OH     +     CH2=CHCH0  (16) 

i-propyl  acraldehyde 
alcohol 

A  n-propoxy  radical  can  react  with  an  allyloxy  radical 
in  two  ways  to  give 

a)  propionaldehyde  and  allyl  alcohol  b)  i-propyl  alcohol 
and  acraldehyde. 
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n-CjHyO-   +  C3H50-   >  CH3CH2CHO     +     CH2=CHCH20H  (17) 

Propion-  allyl  alcohol 

aldehyde 

n-C3H70*  +  C3H50-   >  i-C3H70H     +     CH2=CHCH0  (18) 

i-propyl  acraldehyde 
alcohol 

Two  allyloxy  radical  can  react  to  give  acraldehyde  and 
allyl  alcohol. 

2  C3H50-   >  CH2=CH2CHO     +     CH2=CHCH20H  (19) 

acraldehyde        allyl  alcohol 

Allyl  alcohol  and  acraldehyde  are  formed  by  Reactions 
15,  16,  17,  18  and  19.  The  above  Reactions  11  through  19 
explain  the  formation  of  the  major  products  like  acetone, 
i-propyl  alcohol,  propionaldehyde,  n-propyl  alcohol,  acral- 
dehyde, and  allyl  alcohol.  In  the  reaction  sequence  just 
presented,  conversion  of  two  peroxy  radicals  to  alkoxy 
radicals,  and  subsequent  conversion  of  two  alkoxy  radicals 
to  an  alcohol  and  a  carbonyl  compound  was  shown  as  two 
separate  steps.  A  mechanism  is  known,  however,  which 
combines  these  two  processes.  In  the  "Russel  mechanism", 
two  peroxy  radicals  can  combine  to  give  a  cyclic  transition 
state,  which  subsequently  fragments  to  give  a  carbonyl 
compound,  an  alcohol,   and  molecular  oxgen  (77). 
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2    {R)2CH00-  >  C  0   >    (R)2C0  +  RCH{OH)R 

R^    li  0^ 

I  +  O2 

R— C— R 

I 

H 

Although  this  process  is  interesting  from  the  mechan- 
istic viewpoint,  it  does  not  lead  to  any  new  or  different 
predictions  regarding  product  formation. 

Hot  hydrogen  atoms  born  in  Reactions  2  and  3  may  react 
while  hot,  or  may  become  thermalized. 

H*  +  O2   >  H02-  (20) 

+  C3H8   >  i-C3H7-  +  H2  (21) 

H*  +  C3H8   >  n-C3H7*  +  H2  (22) 

H*  +  C3H8  >  C3H8  +  H-  (23) 

In  most  cases  the  fate  of  thermal  hydrogen  atoms  is 
reaction  either  with  substrate  or  oxygen.  At  this  point, 
however,  the  difference  in  energy  between  primary  and 
secondary  C-H  bonds  becomes  significant;  abstraction  by 
thermal  hydrogen  atoms  gives  predominantly  i-C3H7*. 

H-  +  C3H8   >   i-C3H7-   +  H2  (24) 

H-  +  O2   >  H02-  (25) 

Perhydroxyl  radicals  formed  in  Reactions  20  and  25 
combine  to  give  hydrogen  peroxide,  the  yield  of  which  was 
not  measured. 


99 


H02*  +  H02*  >  H2O2  +  02  -  (26) 

Since  the  yields  of  the  carbonyl  compounds  are  greater 
than  the  corresponding  alcohols,  the  carbonyl  compounds  are 
formed  in  some  other  way  besides  Reactions  11,  12,  13,  14, 
15,  16,  and  19.  A  perhydroxyl  radical  can  also  abstract  a 
hydrogen  atom  from  i-propoxy,  n-propoxy  or  allyloxy  radi- 
cals to  give  acetone,   propionaldehyde  and  acraldehyde. 

i-CsHyO-   +  H02'   >   (CH3)2C=0  +  H2O2  (27) 

acetone 

n-C^UjO'  +  H02'   >  CH3CH2CHO  +  H2O2  (28) 

Propion- 
aldehyde 

C3H50-   +  H02-  >   CH2=CHCH0     +     H2O2  (29) 

Acraldehyde 

Reactions  27,  28,  and  29  explain  why  the  yield  of 
acetone  is  greater  than  iso-propyl  alcohol,  propionaldehyde 
is  greater  than  n-propyl  alcohol  and  acraldehyde  greater 
than  allyl  alcohol. 

All  the  reactions  above  account  for  the  formation  of 
major  products  at  high  concentrations  of  oxygen.  At  low 
concentrations  of  oxygen  it  is  seen  from  Figs.  30a  and  31a 
that  the  yields  of  acetone,  iso-propyl  alcohol  and  n-propyl 
alcohol    increase,    the  yield   of   propionaldehyde  decreases. 
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while  the  yield  of  acraldehyde  and  allyl  alcohol  remains 
the  same.  This  can  be  explained  as  follows: 

At  low  concentrations  of  oxygen  the  hot  hydrogen  atoms 
are  less  efficiently  scavenged  by  oxygen.  Instead,  the  hot 
hydrogen  abstracts  a  hydrogen  atom  from  the  substrate  to 
give  i-C3H7*  and  as  explained  earlier.  The  consequent 

excess  production  of  i-C3H7*  (and  subsequently  i-C3H702* 
and  i-C3H70*)  as  compared  to  n-C3H7*  (plus  n-C3H702*  and  n- 
C3H70*)  accounts  for  enhanced  yield  of  acetone  and  i-propyl 
alcohol  as  oxygen  concentration  falls  off.  The  n-propyl 
alcohol  yield  increases  slightly  at  low  oxygen  concentra- 
tions due  to  the  cross-reaction  with  i-propoxy  radicals. 

It  might  appear  that  the  propionaldehyde  yield  should 
also  increase  at  low  oxygen  concentrations,  due  to  cross- 
reaction  13  and  14.  From  experimental  results,  however,  it 
was  known  that  this  yield  decreases  at  low  oxygen  concen- 
tration. A  large  number  of  computer  simulation  runs  was 
made  in  an  attempt  to  clarify  this  point.  It  was  found  that 
no  reasonable  combination  of  rate  constants  for  reactions 
considered  up  to  this  point  could  explain  the  experimental 
fact  (fall-off  of  propionaldehyde  at  low  oxygen  concentra- 
tion), although  some  combinations  of  parameters  gave  a 
propionaldehyde  yield  which  was  nearly  independent  of 
oxygen  concentration. 
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A  review  of  the  published  literature  on  radiolytic 
oxidation  of  hydrocarbons  provided  a  possible  resolution  of 
the  propionaldehyde  yield  problem.  In  interpreting  the 
results  of  liquid-phase  radiolytic  oxidation  of  n-heptane, 
Kunjappu  and  Rao  (60)  suggested  a  direct  reaction  between 
molecular  oxygen  and  excited  n-heptane  to  give  heptanones. 

C7H16*  +  O2   >  C7H14O  +  H2O 

A  similar  postulate  is  capable  of  explaining  the 
present    results;    a    reaction  of    C3H8*    with  to  give 

propionaldehyde  and  water  is  suggested. 

C3H8*   +  O2   >  CH3CH2CHO  +H2O  (30) 

The  yield  of  this  reaction  would  be  expected  to  fall- 
off  at  low  oxygen  concentration,   since  Reaction  30  must 

compete  with  other  oxygen  independent  channels  for  loss  of 
* 

C3H8  . 

Since  the  concentration  of  allyl  radicals  does  not 
change  at  various  concentrations  of  oxygen,  there  is  not 
much  change  in  the  yields  of  allyl  alcohol  and  acraldehyde 
as  a  function  of  oxygen  concentration.  The  yields  of  the 
latter  products  are  relatively  small  compared  to  acetone, 
propionaldehyde,  iso-propyl  alcohol,  and  n-propyl  alcohol 
in  any  event. 
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Computer  Simulation  of  Radiolytic  Oxidation  Mechanism 

Since  the  entire  reaction  mechanism  too  complex  to  be 
suitable  for  modeling,  several  simplifications  were  made  in 
selecting  the  subset  of  28  reactions  given  in  Table  10  as 
the  basis  for  the  computer  work.  This  scheme  involves  G- 
values  for  only  4  postulated  primary  processes,  along  with 
24  rate  constants  for  subsequent  thermal  and  hot  radical 
reactions.  Unfortunately,  only  8  of  needed  rate  constants 
were  available  from  the  literature  (68,69,71,74-76),  and 
several  of  these  had  a  broad  range  of  uncertainty. 

The  numbers  of  "degrees  of  freedom  "  in  dealing  with 
the  mechanisn  is  considerably  less  than  the  total  number  of 
unknown  rate  constants  involved.  This  is  true  mainly 
because  each  species  has  only  two  or  three  reaction 
channels  by  which  it  is  consumed;  only  the  ratio  of  the 
respective  rate  constants  affects  the  predictions  made  by 
computer  modeling.  There  is  also  some  limitation  on  the 
rate  constants  selected,  in  that  the  numbers  must  be  of 
approximate  magnitude  for  radical-radical  reactions, 
radical-oxygen  reactions,  hydrogen  abstraction  from 
substrate,  etc. 

Several  of  the  simplifications  made  in  choosing  the 
reaction  scheme  for  computer  modeling  were  discussed  above, 
particularly  the  elimination  of  ionic  processes  for  further 
consideration.  Some  additional  simplifications  were  made 
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for  entirely  pragmatic  reasons.  For  example,  cross- 
combination  reactions  between  the  three  peroxy  radicals  (i- 

^3^7^2*'    '^■^3^702*'  ^3^502*)    to   form   the  corresponding 

alkoxy  radicals  and  molecular  oxygen  were  ignored,  since 
including  the  necessary  three  reactions  would  add  little 
new  information.  Additionally,  where  reactions  of  hydro- 
peroxy  radicals  (HOj*)  with  the  three  alkoxy  radicals 
shown,  reactions  of  H02*  with  the  alkyl  peroxy  species  was 
excluded  arbitrarily.  Also,  the  disproport ionation  reaction 
of  2  n-C3H50'  radicals  to  give  the  corresponding  aldehyde 
and  alcohol  is  shown,  but  the  cross-reactions  with  i-C3H70* 
and  n-C3H70*  are  ignored,  as  a  matter  of  convenience.  Even 
with  these  simplifications,  the  entire  scheme  is  rather 
complex. 

More  then  100  computer  runs  were  made  using  various 
combinations  of  mechanistic  reactions  and  their  rate  con- 
stants. The  best  results  were  obtained  by  using  the  combi- 
nation of  mechanistic  steps  and  rate  constants  given  in 
Table  9.  Tables  11  and  12  show  the  experimental  and  compu- 
ter predicted  G-values  of  the  radiolysis  products  at  diff- 
erent oxygen  concentrations.  Table  12  shows  G-values  of 
major  primary  processes,  as  inferred  from  computer  modeling 
work.  Table  13  gives  a  comparison  of  the  measured  and 
calculated  G-values  for  major  products,  as  derived  from 
Tables  11  and  12  for  the  case  of  10%  oxygen.  Figures  28  and 
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Fig.    28     Computer  predicted  yields   of  acetone,  O  ; 

i-propyl  alcohol, B  and  propionaldehyde ,Q 
as   a  function  of  radiolysis  time. 
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29  show  the  computer  predicted  yields  of  major  oxygen 
containing  products  as  a  function  of  radiolysis  time,  over 
a  time  period  of  100  seconds.  The  computations  could  be 
extended  to  a  longer  time,  by  using  more  iterations,  but  it 
was  not  feasible  to  increase  the  time  increment  per  compu- 
ter step,  since  behavior  of  the  concentration  of  transient 
species  become  erratic  if  the  AT  per  computer  increment 
represents  an  excessively  long  time  compared  with  the  half 
lives  of  the  radical  processes.  Figures  30a,  30b,  31a,  and 
32b  compare  experimentally  measured  and  computed  yields  of 
the  oxygenated  products  as  a  function  of  percent  of  added 
oxygen  in  the  radiolytic  oxidation  of  propane.  It  will  be 
seen  that  the  computer  modeling  scheme  provides  a  very 
satisfactory  interpretation  of  the  experimental  data. 

A  typical  printout  made  during  one  run  of  the  numeri- 
cal integration  program  is  given  in  Appendix  V;  it  shows 
the  steady  state  and  net  product,  mechanistic  steps,  their 
rate  constants,  and  the  concentrations  of  several  products 
at  various  time  periods.  Each  print  interval  represents  10 
increments  of  the  integrated  package,  so  that  the  printout 
showing  50  print  intervals  represents  500  iterations  of  the 
program, 

C.  Summary 

The    radiolytic    oxidation    of    gaseous    propane  was 
studied  at  100  torr  pressure  at  various  concentrations  of 
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oxygen  at  25  °C.  The  major  oxygen  containing  products  and 
their  G-values  when  10%  oxygen  was  added  to  propane  were 
acetone,  0.98  ;  i-propyl  alcohol,  0.86;  propionaldehyde, 
0.43;  n-propyl  alcohol,  0.11;  acraldehyde,  0,14;  and  allyl 
alcohol  0.038.  The  minor  products  were  iso-butyl  alcohol, 
tert-amyl  alcohol,  n-butyl  alcohol,  n-amyl  alcohol,  and 
iso-amyl  alcohol.  Small  yields  of  iso-hexyl  alcohol  and  n- 
hexyl  alcohol  were  also  observed.  There  was  no  apperent 
change  in  the  G-values  at  pressures  of  50,  100,  and  150 
torr.  The  formation  of  major  products  was  explained  on  the 
basis  that  the  radicals  H-,  i-C3H7-,  n-C^Hj',  and  CjHs', 
combine  with  to  give  peroxy  radicals  HOj',  i-CjHyOj',  n- 
^3^7°2*'  C3H5O2..    The   peroxy   radicals    i-CjH^Oj,  n- 

^3^702*'  C3H502*  react  with  one  another  to  give  the 

corresponding  alkoxy  radicals  i-C3H70*,  n-C3H70*  and 
C3H50* . 

These  alkoxy  radicals  react  with  one  another  to  give 
carbonyl  compounds  and  alcohol.  Hydroperoxy  radical  also 
reacts  with  the  alkoxy  radicals  to  give  the  corresponding 
carbonyl  compounds.  It  is  suggested  that  propionaldehyde, 
which  is  one  of  the  major  products,  is  also  formed  directly 
by  reaction  of  excited  propane  and  oxygen.  The  reaction 
scheme  for  the  formation  of  major  products,  at  various 
concentrations  of  oxygen,  was  examined  by  computer  modeling 
based  on  a  mechanism  involving  28  reactions.  Yields  could 
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be  brought  into  agreement  with  the  data  within  experimental 
error  in  nearly  all  the  cases. 

One  conclusion  which  emerged  quite  clearly  from  com- 
parative work  on  pure  propane  and  propane-oxygen  systems  is 
the  following:  it  is  not  possible  to  choose  a  single  set  of 
"primary  yields"  which  is  appropriate  for  pure  propane  and 
for  propane  heavily  scavenged  with  oxygen  (1-15%  concentra- 
tion). It  appears  obvious  that  addition  of  oxygen  at  these 
concentration  levels  interferes  with  species  which  act  as 
precursors  to  the  radical  (and  to  a  lesser  extent  molecu- 
lar) yields,  in  the  unscavenged  system.  Nevertheless,  it  is 
of  interest  that  the  total  G-values  for  the  primary  pro- 
cesses are  quite  similar  in  pure  propane  and  propane-oxygen 
system,    7.6  and  8.5  respectively. 

Only  a  single  type  of  propyl  radical  was  assumed  for 
the  pure  propane  system.  It  was  necessary  to  distinguish 
between  normal  propyl  and  iso-propyl  radicals  in  the  case 
of  propane-oxygen  mixtures,  since  the  behavior  of  the  two 
species  is  distinctly  different  when  oxygen  is  added  to  the 
system.  Normal  propyl  radical  is  precursor  to  propion- 
aldehyde  and  n-propyl  alcohol,  while  iso-propyl  radical  is 
a  precursor  to  iso-propyl  alcohol  and  acetone. 

Another  feature  of  this  work  concerns  the  apparent 
absence  of  oxygenated  species  containing  less  than  three 
carbon  atoms.  Similarly  Kunjappu  and  Rao  in  the  radiolytic 
oxidation  of  n-heptane  found  essentially  a  total  lack  of 
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oxygenated  products  containing  fewer  carbon  atoms  than  the 
parents. 


APPENDIX  I 

CALCULATION  OF  DOSIMETRY 

Unlike  a  typical  liquid  phase  radiolysis  where  the 
incident  gamma  rays  interact  directly  with  the  liquid,  the 
usual  case  in  a  gas  phase  radiolysis  is  to  find  that  the 
incident  gamma  rays  interact  mainly  with  the  walls  of  the 
vessel  and  the  ejected  electrons  in  turn  irradiate  the  gas. 
If  the  vessel  meets  certain  geometric  criterion,  which  are 
below,  it  can  be  considered  to  be  a  Brag-Gray  cavity 
(78,79).  This  means  that  the  ionization  in  the  cavity  is 
produced  by  secondary  electrons  and  that  energy  deposition 
in  gases  in  the  cavity  depends  on  their  electron  stopping 
power,    rather  than  their  electron  densities. 

The  nickel  radiolysis  vessels  used  in  the  work  can  be 
shown  to  meet  the  requirements  for  a  Bragg-Gray  cavity  as 
outlined  by  Hine  and  Brownell  (66). 

1.  The  cavity  must  have  dimensions  such  that  only 
a  very  small  fraction  of  the  particle  energy  is 
dissipated  in  it.  This  also  means  that  only  a  very 
small  fraction  of  the  particles  contributing  to 
the  ionization  will  enter  the  cavity  with  a  range 
smaller  than  the  cavity  dimensions. 
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The  radiolysis  vessel  (Fig,  5)  has  an  internal  dia- 
meter of  3.8  cm.  In  a  typical  experiment,  the  ethylene 
pressure  was  300  torr,  which  was  used  as  a  dosimeter.  This 
corresponds  to  a  density  of  4.52x10"^  g/cc.  The  average 
energy  of  compton  electrons  produced  by  Co-60  gamma  rays  is 
0.587  MeV.  These  electrons  have  a  range  of  0.207  g/cm^  in 
ethylene  (80).  Dividing  this  value  by  the  density  of 
ethylene  gives  a  range  of  4.58  meters.  Thus  the  average 
secondary  electron  easily  traverses  the  radiolysis  cell. 

2.  Direct  absorption  of  quantum  radiation  by  the 
gas  in  the  cavity  should  contribute  only  a  negli- 
gible portion  of  the  total  ionization  (66). 

The  gamma  rays  from  Co-60  have  energies  1.17  MeV  and 
1.32  MeV.  In  this  region  of  energies  the  cross  section  for 
the  photoelectric  effect  and  for  pair  production  are  very 
small  compared  to  the  cross  section  for  compton  scattering. 
Thus  the  relative  energy  absorption  of  the  nickel  walls  and 
the  gas  inside  is  proportional  to  the  ratio  of  their  elec- 
tron densities.  The  electron  densities  of  nickel  and 
ethylene  at  300  torr  are  2.56x10^4  electron/cc  and 
1.56x10  electron/cc  respectively.  The  ratio  of  these 
values  is  about  10^.  This  indicates  negligible  absorption 
by  the  gas  compared  to  the  walls  of  the  vessel. 

3.  The  cavity  must  be  surrounded  by  an  equili- 
brium thickness  of  solid  medium  so  that  all  the 
particles  crossing  it  originate  in  the  medium. 
The  required  thickness  for  electronic  equilibrium 
IS    theoretically    equal    to    the    range    of  the 
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particle  of  maximum  energy  released  in  the  medium, 
although  in  practice  path  obliquity  leads  to  the 
approximate  establishment  of  electronic  equili- 
brium at  a  wall  thickness  of  considerably  less 
than  the  maximum  range  (66), 


The  range  of  a  1.3  MeV  electron  in  nickel  is  about 
0.5  mm.  Thus  the  1.19  mm  thick  walls  of  the  vessel  are 
sufficient  to  establish  electronic  equilbrium. 


4.  The  energy  dissipated  by  the  ionizing  particles 
must  be  sensibly  uniform  over  the  volume  of  the 
medium  immediately  surrounding  the  cavity  and  con- 
tributing to  the  electronic  equilibrium  (66). 


This  criterion  should  not  be  critical  in  the  deter- 
mination of  comparative  rather  than  absolute  dose  rates, 
and  all  that  was  necessary  in  the  present  work  was  to  know 
the  total  dose  rate  averaged  over  the  volume  of  the  vessel. 

Now  it  is  possible  to  relate  the  energy  deposited  in 
the  gas  used  as  a  dosimeter  (Ep)  to  that  deposited  in  some 
other  gas  of  interest  (Eg)  by  multiplying  by  the  ratio  of 
their  electron  stopping  powers. 


where  Sg  and  are  the  mass  stopping  powers  of  the  gas 
sample  and  dosimeter  respectively.  The  mass  stopping  power 
is  defined  as 


D 


dE 

S  =  -   

dx 


1 

X  _    erg  cm 

P 
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where  p  is  the  density  of  the  medium  and  -(dE/dx),  usually 
known  as  the  specific  energy  loss  or  stopping  power,  is  the 
energy  loss  per  unit  path  length  by  a  particle  crossing  the 
medium.  This  value  can  be  calculated  by  using  the  Bethe 
equation,   which  has  the  form  (36) 

dE         2lTNe4z  ^oV^e  , 

—  =    {In  _   -   [2(1-32)  ^/^_i+e2]in2 

dx  moV  2l2(l-e2) 

+  1  -  6^  +  i  [1-(1-b2)  ^/2j2  gj.g/^j^ 
8 

When  applied  to  molecules,  the  symbols  have  the  following 
meanings: 

N    =  number  of  molecules/cm^ 
mo  =  rest  mass  of  electron 

e    =  v/c  =  ratio  of  electron  velocity  to  that  of  light 
e    =  charge  on  electron 

Z    =  atomic  number  or  number  of  electrons  per  molecule 

V    =  velocity  of  electron 

I    =  mean  excitation  potential 

Taking  the  ratio  of  the  stopping  power  of  the  sample 
(Ss)  to  that  of  the  dosimeter  (S^)  ^^d  rearranging  the 
terms  in  the  brackets  gives 

Ss  _  (NsZs)  (k-lnl|) 
Sd         (NdZd) (k-lnig) 
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E 


S  ^  Sg  _  (PsZs) (k-lnl§) 
Ed  (PdZd^ (k-lnlg) 


m_  v^E 


«"q  V     u  1/2 

where   k  =  In    -  [2(l-g2)  ^  ^-i+e2)  j  ^^^2  +  1 

2(1-62) 

-   32  +  -  [1.(1-62)'/']' 
8 

Pg  and  Pj)  are  the  pressures  of  the  gas  sample  and 
doeimeter  respectively. 

The  average  energy  of  a  Co-60  gamma  ray  source  is  1.25 
MeV  (average  of  1.17  and  1.33  MeV).  The  average  energy  of 
the  Compton  electrons  generated  by  scattering  of  a  1.25  MeV 
gamma  ray  is 


Eav       —  X  E (photon) 
o 


where  a     is  the  total  Compton  cross  section  per  electron 

Og    is   the   cross   section   for   energy  absorbed  per 
electron 
a    =  1.892xl0"25  cm2/electron 
ag  =  0.899xl0~25  cm^/electron 

(0.899xl0"25) 

Eav  "   ^ —  X  1.25  MeV 

(1.892x1025) 
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=  0.587  MeV 
=  9.404x10"'^  erg 
From  the  relativistic  energy  relation 

Substituting  the  values 
rtiQ  =  9.1095x10-25  g 

c    =  2.9979x1010  cm/sec 
E    =  9.404x10-7  erg 

We  get 

B  =  8.851x10-1 
V  =    c  =  2.653x1010  cm/sec 
Using  the  above  constants,  k  may  be  evaluated 

[m  v^E] 


1  o 
In  

2  (1-  S2) 

-ln2[2(l-B2)l/2_i+g2j  ^  -0.4950 

1-B2 


=  -27.300  (a) 


(b) 

=  0.2166  (c) 


1 


-  [l-{l-62)l/232  ^  0^035^ 


k  =  a  +  b  +  c  +  d 


=  -27.5427 
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Values  for  the  mean  excitation  potentials  (T)  for 
ethylene  and  propane  were  calculated  from  values  for  the 
atoms  (79),   using  the  relation 


(  I  NiZ.  Inl.) 
Inl  =   .  1_ 


where 


=    number  of  atoms  of  "i"  per  molecule 
Zj  =  atomic  number  of  element  "i" 

=  average  excitation  potential  of  element  "i" 
Values  for  1^  can  be  obtained  from  Hine  and  Brownell 
(79)   and  J.  E.  Tumer,  Patricia  D.  Roecklein  and  R.  B.  Vora 
(81)  . 


Element 
Hydrogen 
Carbon 
Oxygen 


Z 
1 
6 
8 


I(eV) 
18.2 
81.3 
101 


I (erg) 

2.91564x10"^^ 
1.30243x10"^^ 
1.61802xl0"10 


Average  values  (T)  for  ethylene,  propane  and  oxygen 
were  55.92  eV,  51.30  eV  and  101  eV,  using  these  values  of 
I,   the  stopping  power  ratio  was  found  by  using  equation 

E(C3H8)       S(C3H8)  P{C3H8) 

  =    =  1.62  X   

E{C2H4)        S(C2H4)  P(C2H4) 
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E(02)  S(02)  P(02) 
  =    =  1.03  X   


E(C2H4)       S(C2H4)  P(C2H4) 

where  P  is  the  pressure  of  the  gas. 

The  total  dose  in  a  mixture  of  gases  can  be  obtained 
by  summing  the  doses  of  the  individual  component 

E(C2H4) 

=  E(C3H8)   +  E{02)   =    x   [1.62P{C3H8)   +  1.03P(O',)] 

P(C2H4) 

To  find  the  G-value  i.e.  number  of  molecules  formed 
per  100  eV  energy  absorbed 


molecules 

G  =    X  100 

eV 


(molecules/gh)   x  100 
G  =  

(eV/gh) 


molecules        G  eV 
gh  100  gh 


moles  G  eV  1 
  =    X           X  — 

gh        100      gh  Nq 

(moles/gh)  x  100  x  N, 


G  = 


(eV/gh) 

where  =  Avogadro's  number  and  G  for  hydrogen  in  ethylene 
dosimetry  is   taken  as  1.2  (64,65). 


APPENDIX  II 


INTERACTION  OF  GAMMA  RAYS  WITH  MATTER 

Gamma  rays  are  electromagnetic  radiations  similar  to 
x-rays,  but  typically  of  shorter  wavelength.  A  significant 
difference  between  x-rays  and  gamma  rays  is  the  source. 
The  source  for  x-rays  is  extranuclear ,  that  of  gamma  rays 
intranuclear,  x-rays  result  from  transitions  between  elec- 
tronic energy  levels,  gamma  rays  on  the  other  hand  are 
associated  with  transitions  between  nuclear  energy  levels. 
Gamma  rays  like  x-rays  are  highly  penetrating,  the  effec- 
tive range  depends  on  the  energy. 

The  absorption  of  gamma  rays  in  matter  occurs  by 
mechanisms  which  are  completely  different  from  the  absor- 
ption of  charged  particles.  The  interaction  is  mainly  with 
the  atomic  electrons,  but  a  gamma  ray  can  lose  a  large 
fraction  of  its  energy  or  all  of  it  in  a  single  encounter. 
The  rays  are  absorbed  according  to  an  exponential  law  which 
is  characterized  by  a  half  thickness  (x^/j)  and  absorption 
coefficient  (p).  if  is  the  initial  intensity  of  a  beam 
of  gamma  rays,  then  the  intensity  (I)  of  rays  of  initial 
energy  after  passing  through    x  cm  of  absorber  will  be 
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The  linear  absorption  coefficient  (v)  is  a  function  of 
the  energy  of  the  gamma  rays  and  the  absorbing  medium. 
Half   thickness    (Xj^/j)  related    to    linear  absorption 

coefficient  by  the  equation 

0.693 

^1/2  "   

V 

Gamma  rays  interact  with  matter  in  several  ways, 
photoelectric  effect,  Compton  effect  and  pair  production 
are  of  major  importance  while  Rayleigh  scattering,  photo- 
disintegration  or  nuclear  reactions  and  nuclear  resonance 
scattering  are  of  less  importance. 

1.  Photoelectric  Effect  :  It  is  the  most  significant 
for  gamma  rays  of  low  energy  and  for  absorbers  of  high 
atomic  number.  The  inelastic  collision  of  the  photon  with 
orbital  electron  results  in  the  complete  ejection  of  the 
electron  and  production  of  an  ion  pair.  The  absorbed 
photon  vanishes.  If  E  is  the  energy  is  the  gamma  ray 
photon,  then  in  a  photoelectric  encounter  an  amount  P, 
equal  to  the  binding  energy  of  the  electron  in  the  atom, 
will  be  required  to  remove  the  electron  and  the  remainder 
(E-P)  is  carried  off  by  the  electron  in  the  form  of  kinetic 
energy 
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K.E.   =  E-P 

(l/2)inv2  =  hv-P 

In  this  process  the  electron  most  probably  comes  from 
the  K  shell,  provided  the  gamma-ray  has  sufficient  energy 
to  supply  the  K  shell  binding  energy.  The  photoelectric 
absorption  follows  the  equation 

Probability  of  photoelectric  absorption   «    z^/(E  ^/^j 

Hence  it  is  most  important  for  absorption  of  low 
energy  gamma  in  heavy  elements. 

2.  Compton  Effect:  This  is  important  for  gamma  rays 
of  medium  energy  (0.5  to  1  MeV)  and  when  the  absorber  is  a 
material  of  low  atomic  number.  It  involves  a  collision 
between  a  photon  and  an  electron  in  which  a  part  of  the 
energy  of  the  photon  is  imparted  to  the  electron.  The 
photon  emerges  from  the  collision  in  a  new  direction  and 
with  reduced  energy.  The  photon  of  lower  energy  may  then 
undergo  a  further  photoelectric  or  Compton  interaction,  or 
the  energy  may  be  so  greatly  diminished  that  it  is  longer 
detectable. 

The  Compton  effect  is  shown  in  Fig.  32.  Considering 
that  both  energy  and  momentum  must  be  conserved  in  the 
coll ision 
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(3) 


(2) 


(1) 


The  first  equation  follows  from  the  conservation  of 
energy.  The  second  equation  represents  the  conservation  of 
momentum  along  the  direction  of  the  incident  photon;  the 
momentum  of  the  electron  is  p,  and  the  momentum  of  a  photon 
of  energy  E  is  E/c.  The  third  equation  expresses  the 
conservation  of  momentum  at  90°  to  the  direction  of  the 
incident  photon.  The  momentum  of  the  electron  is  related 
is  its  kinetic  energy  Eg  by  the  relativistic  in  equation 


where  m  is  the  mass  of  an  electron  at  rest. 

The  above  equation,  with  equations  1,2  and     3  may  be 


~    I  "  (h/mc)  (1  -  costj" ) 

where  and  are  wavelengths  of  the  incident  and 
Compton  scattered  photons.  The  term  h/mc  is  called  compton 
wavelength  and  has  a  value  2.43xl0~10  cm. 

The  mass  absorption  coefficient  due  to  Compton 
scattering  varies  very  little  with  the  atomic  number  of  the 
absorber.  Compton  scattering  makes  only  a  small  contribu- 
tion to  the  absorption  of  high  energy  photons. 


p2    =  (l/c2)    [Eg (Eg  +  2mc2)] 


used  to  eliminate  p.  Eg  and  Q,  with  the  result 
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3.  Pair  Production:  Photons  of  energy  greater  than 
1.02  MeV  may  be  absorbed  by  pair  production.  The  photon 
vanishes  to  be  replaced  by  an  electron  positron  pair  whose 
kinetic  energies  Eg  and  Ep  are  related  to  the  incident 
photon  energy  Ej  by  the  equation 

Ej  =  Eg  +  Ep  +  2mc^ 

The  2mc2  (1.02  MeV)  represents  the  mass  energy  which 
must  be  supplied  to  create  the  particles.  The  process  is 
not  exactly  the  reverse  of  electron-positron  annihilation, 
because  in  that  process  two  photons  are  emitted,  where  as 
in  pair  production,  only  one  is  absorbed.  The  positron  is 
slowed  down  in  a  similar  manner  to  an  electron,  and  even- 
tually combines  with  an  electron,  the  two  particles  being 
replaced  by  two  0.51  MeV  gamma  rays  (annihilation  radia- 
tion)  emitted  in  opposite  directions. 

In  order  to  simultaneously  conserve  momentum  and 
energy  in  pair  production,  a  third  body  (usually  a  nucleus) 
must  be  present.  In  the  annihilation  process,  energy  and 
momentum  are  conserved  simultaneously  by  the  emission  of 
two  photons.  The  importance  of  pair  production  increases 
slowly  for  photon  energies  above  1.02  MeV.  The  probability 
of  pair  formation  increases  with  the  square  of  the  atomic 
number  of  the  absorbing  material.  In  heavy  elements  and 
high  energy  photons,  pair  production  is  much  more  important 
than  the  photoelectric  and  Compton  effects. 
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4.  Rayleigh  Scattering:  If  the  interaction  between  a 
photon  and  an  orbital  electron  is  insufficient  to  produce 
ionization  or  excitation,  the  collision  is  elastic  and  the 
energy  of  the  photon  is  the  same  after  the  collision  as  it 
was  before.  Rayleigh  scattering  rarely  results  in  scatter- 
ing of  the  photon  through  more  than  a  few  degrees. 

5.  Photodisintegration  or  Photonuclear  Reactions: 
Here  the  interaction  involves  a  collision  of  a  high  energy 
photon  (above  8  MeV)  with  a  nucleus.  The  photon  is  comple- 
tely absorbed  in  the  process,  and  a  neutron,  photon  or 
alpha  particle  is  ejected  the  excited  nucleus. 

6.  Nuclear  Resonance  Scattering:  If  the  vibrational 
frequency  of  a  nucleus  is  equal  to  that  of  an  incident 
photon,  obsorption  of  the  photon  may  occur.  The  photon  is 
then  re-emitted  from  the  nucleus.  All  these  interactions 
are  summarized  in  Table  15. 
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Table  15 

Interaction  of  Radiation  with  Matter 


Incident  In  Collision  Type  of  Collision 
Radiation  with   


Elastic      Inelastic  Complete 

Absorption 


Gamma  Nucleus  Nuclear  Photo- 

Resonance  disintegration 


Orbital  Raleigh  Compton  Photoelectric 
Electron    Scattering  Effect  Effect 


Field  Pair 

Production 
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Incident  Photon,  E, 


Compton  Scattered 
Photon,  E 


Ejected  Electron,  E 


Electron 


Incident  Photon 
 ^  


Pos  itron 


Annihilation 
Radiation 


Fig.    32     The  Compton  effect  and  Pair-production 


APPENDIX  III 


RELATIVE  FLAME  IONIZATION  DETECTOR  RESPONSE  OF 
PRODUCTS  IN  THE  RADIOLYSIS  OF  PROPANE  AND  IN 
THE  RADIOLYTIC  OXIDATION  OF  PROPANE 

In  the  radiolysis  of  propane,  prior  to  each  analysis, 
the  response  of  the  flame  ionization  detector  was  deter- 
mined using  gaseous  propane  as  standard.  Flame  ionization 
detector  is  not  sensitive  to  hydrogen;  it  was  measured 
using  thermal  conductivity  detector. 

In  the  radiolytic  oxidation  of  propane,  acetone  was 
used  as  a  standard.  The  relative  flame  ionization  detector 
response  on  Tracer  550  gas  chroma tograph  are  given  in 
Tables  16  and  17. 
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Table  16 

Flame  Ionization  Detector  Response  of 
Products  Relative  to  Propane  on 
Tracor  550  Gas  Chromatograph 


Methane 

V)  »i  1  z 

Ethane 

a  "Til 

Propane 

i.  •aV 

I-Butane 

1.24 

N-Butane 

1.24 

I-Pentane 

1.71 

N-Pentane 

1.71 

Hexane 

2.13 

Heptane 

2.24 

Octane 

2.41 

Nonane 

2.52 

Decane 

2.67 
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Table  17 


Flame  Ionization  Detector  Response  of 
Products  Relative  to  Acetone  on 
Tracor  550  Gas  Chromatograph 


Propionaldehyde 

0.73 

Acetone 

1.00 

Iso-Propyl  Alcohol 

1.39 

N-Propyl  Alcohol 

1.39 

Acraldehyde 

1.13 

Allyl  Alcohol 

1.38 

Iso-Butyl  Alcohol 

2.41 

N-Butyl  Alcohol 

2.75 

Tert-Amyl  Alcohol 

2.69 

Iso-Amyl  Alcohol 

3.69 

N-Amyl  Alcohol 

3.60 

APPENDIX  IV 


DETAE  KINETICS  PROGRAM 


DIMENSION  AK(34) ,C(30) ,CPDHAM(30) ,D(54) ,DATE(2) , 

healerCi 0) 

DIMENSION  ISSd  6),JRE(34,4)  ,JSS(1 6,34)  ,JTENTY(8), 
TBLCPD(8) 

DIMENSION  CH0LD(8, 50) ,TIME(50) 
DIMENSION  TBLH(8) ,DL(34) ,EQUAT(6) ,PE0BLM(12) 
MXAPP=33 
MXCPDS  =  30 
IQ79=0 
MXR  =  34 
MXSTST  =  16 
MXTBNT  =  8 
READ   (7,7000)  BL 
7000  FORMAT  (A4) 

READ   (7,8009)  DAMP.SCALEF 

8009  FORMAT  (2E8.4) 

5  DO  1 0  N  =  1 ,  MXR 
DO  1 0  J  =  1 ,  4 
10  JRE    (N,J)    =  0 

1 2  DO   1 5   J   =  1 ,  MXSTST 

13  ISS    (J)   =  0 

1 4  DO   1 5   I   =  1 ,MXAPP 

1 5  JSS    (J, I)  =  0 

DO  1 7  J   =  1 ,  MXTBNT 
17  TBLH   (j)   =  BL 
20  READ   (7,8000)  DATE, HEADER 
8000  FORMAT  (lX,12A4) 

25   WRITE    (5,9000)   DATE, HEADER 

9000  FORMAT  (12A4) 

30  READ    (7,8010)   NCPDS,    NSS,    NR ,    NRCAT,  NCPCAT 
WRITE  (5,6995) 

6995   FORMAT  (1H2) 

WRITE    (5,8010)   NCPDS, NSS, NR, NRCAT, NCPCAT 

8010  FORMAT  (l4(l4)) 

31  IF    (NCPDS-MXCPDS)  35,35,32 

32  NCPDS   =  MXCPDS 

33  WRITE    (5, 9001  ) 

9001  FORMAT    ('  '  TOO  MANY  COMPOUNDS') 
35  IF  (NSS-MXSTST)  39,39,36 
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56  NSS   =  MXSTST 
37  WRITE  (5,9002) 

9002  FORMAT   ('         '  TOO  MANY   SS  COMPOUNDS') 

39  IF  (NR-MXR)  45,45,40 

40  NR  =  MXR 

41  WRITE  (5,9003) 

9003  FORMAT   C         '  TOO  MANY  REACTIONS') 

45   READ   (7,8020)    ( CPDNAM   (k),   K  =  1,  NCPDS) 
8020  FORMAT   (7 ( 4X , A4 ) , / , 7 ( 4X , A4  )  ) 

WRITE    (5,8020)    (CPDNAM(K) ,K=1 , NCPDS) 
50  DO  77  N=1 ,NR 

52  READ    (7,8010)    (JRE(N,I),  1=1,4) 

53  IF   (NSS)  77,77,54 

54  DO  75  I  =  1 ,4 

55  IF  (jre  (n,i))  75,75,56 

56  IF   (JRE   (N,l)  -  NSS)  58,58,75 
58  J   =   lABS    (jre  (N,I)) 

60  IS   =  ISS    (j)  +1 

61  ISS    (j)   =  IS 

63  IF   (IS  -  MXAPP)  67,67,64 

64  ISS    (J)   =   ISS   (J)   -  1 

65  WRITE   (5,9005)   CPDNAM  (j) 

9005   FORMAT    ('   ',    '  TOO  MANY  APPEARANCES   OF   ',    A4 ) 

66  IS  =  MXAPP 

67  JSS    (J,IS)   =  N 

68  JRE    (N,I)   =   -JRE  (N,I) 
70  IF  (1-3)  72,75,75 

72  JSS    (J,IS)   =  -N 
75  CONTINUE 
77  CONTINUE 
80  NSSI   =  NSS  +  1 
85   IF   (NSS)  110,110,100 
100  WRITE    (5,9010)    (CPDNAM(k),K   =  1,NSS) 
9010   FORMAT    ('   ',   3X,    '   STEADY   STATE   INTERMEDIATES- ', 1 X , 
17  (A4,2X)) 

110  WRITE   (5,9020)    (CPDNAM(k),   K  =  NSSI,  NCPDS) 
9020   FORMAT    ('   ','REACTANTS   AND  PRODUCTS-' ,5X , 7 ( A4 , 2X ) /31 X . 

17(A4,2X)) 
1 20     DO   200  N   =   1 ,NR 

JJ   =   lABS    (JRE    (N,1  )  ) 

EQUATd  )   =   CPDNAM  (jj) 
130   IF   (JRE    (N,2))  140,150,140 
140  JJ   =  lABS    (JRE  (N,2)) 

EQUAT(3)    =   CPDNAM  (jj) 

GO  TO  160 
1  50  EQUAT(2  )   =  BL 

EQUAT(3)   =  BL 
1  60  JJ   =   lABS    (JRE(N,3) ) 

EQUAT(4)   =   CPDNAM  (jj) 

IF   (JRE   (N,4))  170,180,170 
170  JJ   =  lABS    (JRE  (N,4)) 
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EQUAT(6)   =  CPDHAM  (jj) 

GO  TO  182 
180  EQUAT(5)   =  BL 

EQUAT(6)   =  BL 
182   IF   (JRE   (N,2)    .EQ.   O)   GO  TO  184 

IF   (JRE   (N,4)    .EQ.   0)   GO  TO  188 

WRITE  (5,9050)  N,EQUAT(1),  EQUAT(3),  EQUAT(4),  EQUAT(6) 
GO  TO  200 

9030  FORMAT    ('         13,   IX,   A4 ,    '  +         A4 ,   2X ,    '=',2X,   A4 , 
1'   +  A4) 

184  IF   (JRE   (N,4)    .EQ.   O)   GO  TO  186 

WRITE   (5,9031)   N,   EQUAT(I),   EQUAT(3),  EQUAT(4), 
1 EQUAT(6) 
GO  TO  200 

9031  FORMAT   ('         13,   1X,   A4,   3X,   A4 ,   2X ,    '=',   2X ,   A4 , 
1'  +  A4) 

186  WRITE   (5,9032)  H,   EQUAT(I),   EQUAT(3),  EQUAT(4), 
1 EQUAT(6) 
GO  TO  200 

9032  FORMAT    ('         13,   IX,   A4 ,   3X ,   A4 ,   2X ,    '=',   2X ,   A4 ,   3X , 
1  A4) 

188  WRITE   (5,9033)  N,   EQUAT(i),   EQUAT(3),   EQUAT(4),  EQUAT 
1  (6) 

9033  FORMAT    ('         13,    IX,   A4,    '  +         A4 ,   2X,    '=',2X,   A4 , 
13X,  A4) 

200  CONTINUE 

230  IF   (NCPCAT  *  NRCAT)  240,240,235 
235  WRITE   (5,9045)  NRCAT,    CPDNAM  (NCPCAT) 
9045   FORMAT   ('         'REACTION         12,   IX,    'IS  CATALYZED  BY', 
1  A4) 

240  READ   (7,8010)  NPROB,NREPSS 

260  DO  900  NPRB  =  1 ,  NPROB 

270  READ   (7,8000)  PROBLM 

280  READ   (7,8050)    (AK   (n),   N   =  1,NR) 
8050  FORMAT  (7E10.5) 

285  WRITE   (5,9000)   HEADER,  DATE 
WRITE  (5,6995) 

290  WRITE    (5,9060)  PROBLM 
9060  FORMAT    ('   ',    '   ',  12A4) 

300  READ   (7,8060)   NP , NI , DT , T , IWRITE , NTENT 

301  WRITE   (5,8060)  NP , NI , DT , T , IWRITE , NTENT 
8060  FORMAT   ( 15 , 14 , 2E1 0 . 5 , 13 , 14  ) 

305  IF   (NTENT-MXTBNT)  307,307,306 

306  NTENT  =  MXTBNT 

307  READ   (7,8010)    (JTENTY   (l),    1=1,  NTENT) 
WRITE  (5,6995) 

310  WRITE   (5,9070)     DT,    (N,   AK   (n),   N  =  1 ,   NR ) 

9070   FORMAT    ('  DT   =   ' , 1 PE1 3 . 5 , / ,  '   RATE   CONSTANTS-' , 1 X , 

1 5(I2,1X,E1  1  .5,3X)/17X,5(I2,1X,E1 1 .5,3X) ,/,17X, 

25(I2,1X,E1 1 .5,3X),/,17X,5(I2,1X,E11 .5,3X),/,17X, 

35  I2,1X,E1  1  .5,3X),/,17X,5(I2,1X,E1 1  .5,3X),/,17X, 
45(12, IX, El  1  .5,3X)) 
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315  READ   (7,8050)   (c(l),   1=1,  NCPDS) 
READ(7,8050)  CSSTOT 

317  DO  319  I   =  1 ,  HTENT 

318  J   =  JTEKTY  (l) 

319  TBLH   (I)   =  CPDNAM  (j) 
WRITE  (5,6996) 

6996  FORMAT   (1  HI  ) 

320  WRITE    (5,9080)      (TBLH    (i),    1=1,  NTENT) 
9080  FORMAT    ('   ',/,   6X ,    'TIME',    1  OX ,   8(A4,    9X ) ) 

322  IF   (HSS)  332,332,323 

323  DO  326  I  =  1 ,  NSS 

324  IF  (C(l))  325,325,326 

325  C(l)   =  1 .OE-20 

326  CONTINUE 

332  DO  334  1=1,  NTENT 

333  J   =   JTENTY  (l) 

334  TBLCPD    (I)   =  C(J) 

340  WRITE  (5,9090)  T,  (TBLCPD  (i),  1=1,  NTENT) 
9090  FORMAT   ('   ',  1P9E13.4) 

DO  350  N  =  1 ,  NR 
345  DL   (N)   =  0 
350  AK   (n)   =  AK    (h)   *  DT 
360  DO  850  IP  =  1 ,  NP 
370  DO  800  IT  =  1 ,  NI 
375  IF  (nss)  580,580,380 
380  DO   550  IREP   =  1  ,NREPSS 
390  DO  420  N  =  1  ,  NR 

N1    =   lABS    (JRE    (N,1  )  ) 

N2   =   lABS    (JRE  (N,2)) 

IF   (N2)  400,400,410 
400  d(n)  =  ak(n)  *  C(N1 ) 

GO  TO  420 
41 0  D(N)   =  AK(N)  *  C(N1  )   *  C(N2) 
420  CONTINUE 

ITEST  =  0 
430  DO  530  IS  =  1 ,  NSS 

1ST  =  ISS(IS) 

XNUM  =  0. 

DEN  =  0. 
440  DO  470  JS  =  1 ,  1ST 

445  KS  =  JSS    (IS,  JS) 

446  IF  (KS)  460,470,450 
450  XNUM  =  XNUM  +  D(KS) 
452   GO  TO  470 

460  KS  =  -KS 

465   DEN  =  DEN  +  D(KS) 

470  CONTINUE 

IF   (DEN)  472,475,472 
472   F  =  XNUM/DEN 

IF   (F)  477,475,477 
475  F  =  1 . 
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GO  TO  488 
477  IF   (F  -  1 . )480,488,484 
480  F  =  DAMP/F  +  1.-DAMP 

F  =  1 ./F 

GO  TO  488 
484  F  =  F*DAMP  +  1 .  -DAMP 
488  C(IS)   =  C(IS)  *  F 

FV  =  F  -  1  . 
490  IF   (I¥RITE)51 0,51 0,500 

500  WRITE    (5,91  00)lP,IT,IREP,XIIUM,DEN,F¥, 
1 (CPDNAM(l),C(l),I=1 ,NSS) 
9100   FORMAT    ('        'PRINT   INTERVAL' , 13 , 2X , 'ITERATION' , 13 , 2X , 
1  'SS  REPEAT ', 13, 2X, 'NUMERATOR ',1PE1 2. 5, 2X, 
2'DEN0MINAT0R'E12.5,2X,'PER  CENT   CORRN . ' , 2PF1 0 . 3 , / ,  1  X , 
36HCONCS. ,1X,A4,1PE1 2.5,2X,A4,E1 2.5,2X,A4,E1 2.5,2X,A4, 
4E1 2.5,2X,A4,E1 2.5, /A4,1  PE1 2 .5,2X,A4,E1 2.5,2X,A4,E1  2 .5, 
52X,A4,E1 2.5) 

510  IF  (abs(fv)  -  .005)530,530,520 

520  ITEST  =  1 

530  CONTINUE 

531  IF   (CSST0T)540, 540, 532 

532  CS  =  0. 

BO  533   I   =  1 ,  NSS 

533  CS  =   CS  +  C(l) 
F  =  CSSTOT/CS 

DO  537  I   =  1 ,  NSS 
537   C(l)   =  C(l)  *  F 
540  IF   (ITEST)  550,560,550 
550  CONTINUE 

553  WRITE  (5,9110)  NREPSS 
9110   FORMAT    ('   ', 'STEADY   STATE   CALCULATION  DID  NOT 
1  CONVERGE','   IN', 12, 1X, 'ITERATIONS') 
555  WRITE   (5,9100)   IP , IT , IREP ,XNUM,DEN , FW, 
1  (CPDNAM(I) ,C(l) ,1  =  1  ,NCPDS) 
GO  TO  900 
560  IF   (IP  *  IT  -  1 )  580,570,580 
570  DO  575  I   =  1  ,  NTENT 
572  J   =  JTENTY(l) 

575  TBLCPD(I)   =  C(J) 

576  WRITE   (5,9090)   T, (TBLCPD(i),I   =  1,NTENT) 
580  CONTINUE 

585  DO  620  N  =  I ,  NR 

N1    =   IABS(JRE(N,1  )) 

N2  =  IABS(JRE(N,2)) 
590  IF  (N2)  600,600,610 
600  d(h)  =  ak(n)  »  c(N1 ) 

GO  TO  620 
610  D(n)   =  AK(n)   *   C(N1  )   *  C(N2) 
620  CONTINUE 

625  IF  (NCPCAT  *  NRCAT)  630,630,627 
627  D(NRCAT)  =  D(NRCAT)  *  C(NCPCAT) 
630  DO  638  N  =  1 ,  NR 


136 


651   IF  (DL(n))  636,636,632 

632  BLAST  =  D(n) 

633  D(N)  =  1.5  *  D(n)  -  0.5  *  DL(n) 

634  DL(n)   =  DLAST 

635  GO  TO  638 

636  DL(n)   =  D(n) 
638  CONTINUE 

640  DO  680  N  =  1 ,  NR 
642  DO  680  1=1,4 

644  IF   (JRE(N,I))  680,680,645 

645  II   =  JEE(N,I) 

650  IF  (I  -  3)  655,665,665 
655  C(II)  -  C(II)  -  d(n) 
660  GO  TO  680 
665  C(II)   =  C(II)  +  D(n) 

680  CONTINUE 

681  NGCHK  =  0 

682  DO  685  1=1,  NCPDS 

683  IF   (C(l))  684,685,685 

684  NGCHK  =  NGCHK  +  1 

685  CONTINUE 

690  IF   (NGCHK)  750,750,700 
700  WRITE   (5,9120)   IT , ( C ( I ) , I  =  1  , NCPDS ) 
9120  FORMAT    ('  NEGATIVE   CONCS .-( INTERVAL' , 13 , 1 X , ' 

1 C-ARRAY',1P6E14.5,/,40X,6E1  4.5) 
710  WRITE   (5,9130)    (D(I),   I   =  1,  NR) 

9130  FORMAT   ('   ',33X,'D-ARRAY',1PE1 4.5,/,40X,6E14.5) 

720  GO  TO  900 
750  T  =  T  +  DT 
800  CONTINUE 

805  DO  807  1=1,  NTENT 

806  J   =  JTENTY(l) 

807  TBLCPD(I)   =  C(J)*SCALEF 

810  WRITE   (5,9090)   T, (TBLCPD(i) ,1-1 , NTENT) 
81 5  TIME(IP)   =  T 
820  DO  830  KK   =   1 , NTENT 
830  CH0LD(KK,IP)   =  TBLCPD(KK) 
850  CONTINUE 
9095  FORMAT   ('  1P8E13.4) 
860  WRITE   (8,9095)  CHOLD 
900  CONTINUE 
END 


APPENDIX  V 


OUTPUT"  OF  COMPUTER  SIMULATION  PROGRAM 


Fortran  Data  File 


0.7,1 .0 

PROPANE  RADIOLYSIS 

20,9,25,0,0 

H«  H.  CH3.        ET*         C2H5        C3H7  PE* 

C3H5        CTJS*       H2  CH4         C2H6        C4HX  CSHT 

C6HZ       C2H4        C3H6        W  POLY  C5H8 

20,0,9.0 

20,0,6,1 

20,0,10,17 

20,0,11  ,16 

20,0,4,3 

20,0,7,2 

7,0,8,10 

9,0,3,5 

4,20,12,6 

1  ,20,10,6 

1 ,20,2,20 

2,20,10,6 

3,20,11,6 

5.20,12,6 

3.3,12,0 

5.5.12,16 

5.5.13,0 

3,6,13,0 

5,6,14,0 

6,6,15.0 

6,6,20,17 

8,6,17,17 

16,18,19,0 

17,18,19,0 

20,18,19,10 

1  ,400 
PROBLEM  ONE 

3 . 9E-9 , 4 .35E-9 . 6 . 52E-9 , 3 .24E-9 , 4 . 76E-9 , 7 .8E-1  0 , 1  . OEl  2 
9.0E7,1E12,1  .1E9,1  .0E9.9.0E4,1 ,0E2,1  .0E2 
2.5E10,1E10,1  .4E10,1  .2E10,1  .60E1 0,3.4E9,2 .00E9 
2E9,5.0E3,1  .0E4,7.50E-4 
50,20,1E1  ,0.0,0,8 
10,11  ,12,13.14.15.16,17 

1  .OE-20,1  .OE-20,1  .OE-20,1  .OE-20,1  .OE-20,1  .0E-20,1E-20 
1 .OE-20,1 .OE-20.1 .OE-20,1 .OE-20,1  .OE-20,1 .0E-20,1E-20 
1E-20,1E-20,1E-20,1  ,1  E-20,5.38E-3 
0.0,0,0,0 
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Output   of  Computer  Simulation  Program; 
Mechanism  and  System  Parameters 


PEOPARE  RADIOLYSIS 

20     9    25      0  0 

H*           H.  CK3.        ET*         C2K5  C3H7 

C3H5       C3e»  H2          CH4         C2H6  C4HX 

C6HZ       C2H4  C3H6        ¥            POLY  C3H8 

STEADY  STATE  INTEHMEDIATES-  H*       H.        CH3 . 

STEADY  STATE  INTERMEDIATES 


REACTANTS  AND  PRODUCTS- 


C3H5  C38» 
H2       CH4      C2H6    C4ffl  C5HY 
C3H6    W        POLY  C3H8 


PR« 

C5HY 

ET»      C2H5    C3H7  PE* 


C6HZ  C2K4 


1 

C3H8 

•  C38* 

2 

C3H6 

•  C3K7 

H* 

3 

C3He 

-  H2 

C3H6 

4 

C3H8 

•  CH4 

C2H4 

5 

C3H8 

-  ET* 

CH3. 

6 

C3He 

-  PP* 

♦ 

H. 

7 

PR* 

■  C3H5 

K2 

8 

C3e» 

■  CH3. 

* 

C2H5 

9 

ET* 

C3H8 

■  C2H6 

4- 

C3H7 

10 

H* 

C3H8 

■  K2 

♦ 

C3H7 

1 1 

H* 

C3H8 

■  H. 

C3He 

12 

H. 

C3H8 

'  H2 

+ 

C3H7 

13 

CK3. 

C3H8 

•  CH4 

C3H7 

14 

C2H5 

♦ 

C3H8 

■  C2H6 

C3H7 

15 

CH3. 

+ 

CH3. 

■  C2H6 

16 

C2H5 

C2H5  ■ 

■  C2H6 

C2H4 

17 

C2H5 

C2H5  ' 

'  C4HX 

IB 

CH3. 

C3H7  ■ 

■  C4HX 

19 

C2H5 

C3K7  ■ 

■  C5HY 

20 

C3H7 

+■ 

C3H7  ■ 

■  C6HZ 

21 

C3H7 

C3H7  ■ 

■  C3H8 

C3H6 

22 

C3H5 

C3H7  ' 

■  C3H6 

♦ 

C3H6 

23 

C2H4 

+ 

¥ 

•  POLY 

24 

C3H6 

W 

POLY 

25 

C3H8 

W 

POLY 

H2 

PROPANE  RADIOLYSIS 

PROBLEM  ONE 
50  20.1O00OE+02.0OO00E*01 


DT  -      1  .0OO0OE*01 
RATE  CONSTANTS-  1 
6 
11 
16 
21 


3.90000E-09 
7.eOOOOE-10 
1  .000O0E+O9 
1  .OOOOOE*10 
2.00000E*09 


22 


2  4.35000E-09 
7  1  .OOOOOE+12 
12  9.00O0OE*O4 
17  1  .40000E+10 
2.00000E*09 


3  6.520OOE-09 
8  9.000O0E*07 
13  1  .00000E*02 
18  1  .20000E*10 
23  5.0OOO0E+O3 


4  3.24000E-09 
9  1 .OOOOOE+12 
14  1 .0OOO0E*02 
19  1  .60000E*10 
24  1 .0O00OE*O4 


5  4.76000E-09 
10  1.1  OOOOE+09 
15  2.50000E+10 
20  3.40000E+09 
25  7.50000E-04 
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Fortran  Data  File 


0.7,1  .0 

PEOPANE  -  OXYGEN  SYSTEM 

29,10,32,0,0 


H* 

H. 

PR* 

AL. 

IP. 

N?. 

C38* 

IPOD 

NPOD 

ALOD 

H02. 

IP02 

NP02 

AL02 

IPO. 

NPO. 

ALO. 

H2 

H202 

C12P 

H20 

PK 

IPOH 

PA 

NPOH 

ALA 

ALOH 

02 

C3H8 
29,0,7,0 
29,0,5.1 
29.0,6.1 
29,0.3,2 
3.0.4.18 
7.0,20.0 
7,28,24,21 
7,28,29,28 
1 ,28.11 .0 
1,29,18,5 

I  ,29,18,6 
1,29.2,29 
2,29,18,5 
5,28,12,0 
6,28,13,0 
4,28,14,0 
2,28,11  ,0 
12,12,8,28 
8,0,15,15 
13,13,9,28 
9,0,16,16 
14,14,10,28 
10,0,17,17 

II  ,11  ,19,28 
15,15.22,23 
15,16,22.25 
15.16,24,23 
16,16,24,25 
17,17,26,27 
15,11  .22.19 
16,11  .24.19 
17,11  ,26.19 
1  ,400 

PROBLEM  ONE 

4.54E-9,6.40E-9,1 .93E-9,8.86E-1 0,1 .0OE12,9.0OE7,3.25E10 
3.00E10,1  .50E10,1  .00E9,1  .00E8,1  .00E9,8,20E4,4.63E1 0 
2.00E10,1  .0OE10,3.00E1 1  ,1  .OOEIO.I  .00E12,1  .0OE9,1  .00E12 
6.00E10,1  .0OE12,2.0OE9,3.0OE10,1  .00E9,1  .00E9.1  .0OE8 
1 .2E10,9.4E8,6.4E8,4.0E9 
50,10,0.200,0.0,0,8 
18,21  ,22,23,24,25,26,27 

1  .OE-16,1  .OE-16,1  .OE-16,1  .0E-16,1  .0E-16,1  .0E-16,1  .OE-162 
1  .OE-16,1  .OE-16,1  .OE-16,1  .OE-16,1  .OE-16,1  .OE-16,1  .OE-16 
1  .OE-16,1  .OE-16,1  .OE-16,1  .OE-16,1  .OE-16,1  .OE-16,1  .OE-16 
1  .OE-16,1  .OE-16,1  .OE-16,1  .OE-16,1  .OE-16,1  .OE-1 6, 5.38E-4 
5.38E-3 
0,0,0,0,0 


141 


Output  of  Computer  Simulation  Program: 
Mechanism  and  System  Parameters 


PEOPABE  -  OXYGEN  SYSTEM 


29    10    32  0 

0 

H* 

H. 

PE* 

AL. 

IP. 

NP. 

C38« 

IPOE 

NPOD 

ALOD 

H02. 

IP02 

irP02 

AL02 

IPO. 

NPO. 

ALO. 

H2 

H202 

C12P 

K20 

PK 

IPOK 

PA 

NPOK 

ALA 

ALOH 

02 

C3He 

STEADY 

STATE 

INTEEKEDIATES- 

H»  H. 

PE» 

AL. 

STEADY 

STATE 

INTERMEDIATES- 

IPOD  IIPOD 

ALOD 

EEACTABTS  AND  PEODDCTS-         H02 .     IP02    NP02    AL02     IPO.    NPO.  ALO. 

H2        H202    C12P    H20      PK        IPOH  PA 
REACTANTS  AND  PSODUCTS-         NPOH    ALA      ALOH    02  C3H8 


1 

C3HB 

-  C38* 

2 

C3H8 

-  IP. 

H* 

3 

C3H8 

-  NP. 

♦ 

H* 

4 

C3H8 

•  PE* 

H. 

c 

5 

PE« 

-  AL. 

* 

K2 

6 

C7e» 

-  C12P 

7 

C38* 

02 

•  PA 

K20 

8 

C3e« 

+ 

02 

-  C3H8 

02 

9 

w* 

02 

-  H02. 

10 

H» 

+ 

C3H8 

-  H2 

IP. 

1 1 

H* 

♦ 

C3H8 

-  K2 

+ 

NP. 

12 

H* 

C7H8 

-  H. 

C3H8 

13 

H. 

+ 

C3H8 

-  H2 

IP. 

U 

IP. 

02 

•  IP02 

15 

NP. 

02 

•  NP02 

16 

AL. 

02 

-  AL02 

17 

H. 

02 

•  H02. 

18 

IP02 

IP02 

•  IPOD 

02 

19 

IPOD 

•  IPO. 

IPO. 

20 

NP02 

NP02 

-  NPOD 

02 

21 

NPOD 

•  NPO. 

NPO. 

22 

AL02 

AL02 

■  ALOD 

+ 

02 

23 

ALOD 

-  ALO. 

4 

ALO. 

24 

H02. 

H02. 

-  H202 

♦ 

02 

25 

IPO. 

* 

IPO. 

'  PK 

IPOH 

26 

IPO. 

NPO. 

-  PK 

NPOH 

27 

IPO. 

+ 

NPO. 

«  PA 

IPOH 

28 

NPO. 

NPO. 

-  PA 

NPOH 

29 

ALO. 

ALO. 

-  ALA 

ALOH 

30 

IPO. 

H02. 

-  PK 

H202 

31 

NPO. 

H02. 

•  PA 

H202 

32 

ALO. 

K02. 

*  ALA 

H202 

PROPANE  - 

OXYGEN  SYSTEM 

PEOBLEK  ONE 

50    10.20O0OE*00.0OOO0E*01    0  8 

DT  ■  2.00000E-01 

RATE  CONSTANTS-    1  4.540OOE-09  2 

6  9.00000E*07  7 

11  1 .0O000E*O8  12 

16  1  .0OO00E*10  17 

21   1.00000E*12  22 

26  1 .00000E*09  27 

31   6.40000E+08  32 


6.40000E-09 
3.25OO0E*10 
1  .OOOOOE+09 
3.00000E*1 1 
6.00000E+10 
1 .OOOOOE+09 
4.000O0E+09 


3  1 .93000E-09 
8  3.00000E+10 
13  8.20000E+04 
16  1  .OO00OE*10 
23  1 .0O0OOE*12 
28  1 .OOOOOE+08 


4  8.86000E-10 
9  1 .50000E+10 
14  4.630OOE*10 
19  1  .00000E*12 
24  2.000OOE*O9 
29  1 .20000E*10 


5  1 .OOOOOE+12 
10  1 .OOOOOE+09 
15  2.00000E*10 
20  1 .OOOOOE+09 
25  3.00000E+10 
30  9.40O0OE+Oe 
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